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ABSTRACT 


Langevin  equations  for  intrinsic  optical  bistability  in  a  system  of  two-level 
atoms  interacting  by  dipole-dipole  interactions  are  treated,  and  the  stochas¬ 
tic  effects  of  quantum  noise  are  analyzed.  For  cases  in  which  the  population 
relaxation  is  very  slow  (relative  to  the  homogeneous  dephasing  rate  of  the 
complex  dipole)  we  derive  a  one-dimensional  Ito  stochastic  differential  equa¬ 
tion  for  the  inversion  of  population.  For  such  cases  the  probability  distribu¬ 
tion,  diffusion  process  and  the  first  passage  times  between  the  two  metastable 
states  are  calculated.  On  the  other  hand,  the  quantum  fluctuations  in  in¬ 
trinsic  bistability  of  a  pure  quantum  mechanical  system  are  calculated  by 
using  the  Einstein  relations.  We  find  that  a  quantum  noise  which  is  above 
the  standard  quantum  limit  is  inherent  in  the  reaction  field.  We  develop  a 
theoretical  analysis  for  the  spectrum  of  the  light  scattered  by  two,  two-level 
atoms  driven  near  resonance  by  an  external  electromagnetic  field. 


Keywords: 

Intrinsic  Optical  Bistability;  Cooperative  Effects;  Resonance  Fluorescence; 
Quantum  Noise. 
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1.  INTRODUCTION 


Intrinsic  optical  bistability  (10B)  that  is  not  caused  by  external  feedback 
such  as  mirrors,  has  been  the  subject  of  recent  intense  interest.1  It  was  first 
pointed  out  by  Bowden  and  Sung2"  and  by  Bowden26  that  optical  bistability 
(OB)  may  occur  for  a  system  comprising  a  collection  of  atoms  which  inter¬ 
act  with  the  electromagnetic  field  and  are  driven  by  an  externally  applied 
coherent  field  without  external  feedback.  The  first  detailed  experimental 
study  of  intrinsic  optical  bistability  was  conducted  almost  simultaneously  by 
Hajto  and  Janossy3  using  amorphous  GeSe2 ,  who  interpreted  their  results 
as  due  to  temperature-dependent  induced  optical  absorption  in  the  mate¬ 
rial,  and  Bohnert,  Kalt  and  Klingshirn4,  who  used  CdS,  and  also  Rossmann, 
Hennbberger  and  Voigt5  using  the  same  material.  The  process  in  the  first 
case  depends  upon  absorption  due  to  temperature  variation  induced  by  the 
incident  field,  whereas  the  later  cases  depend  upon  saturation  of  absorption 
due  to  the  generation  of  carriers  in  the  material,  and  the  IOB  has  been 
interpreted8  as  due  to  bandgap  renormalization.  Since  the  earlier  works, 
there  have  been  many  theoretical  and  experimental  investigations  of  various 
forms  of  IOB1. 

In  certain  recent  works7  the  nonlinear  oscillator  model  has  been  used  as 
a  prototype  medium  exhibiting  intrinsic  bistability.  Various  effects  related 
to  intrinsic  bistability  have  been  explored  by  following  this  model.  On  the 
other  hand  the  two-level  system  has  been  used  as  a  prototype  medium  for 
describing  quantum  mechanical  effects  in  intrinsic  optical  bistability.8-12  In 
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previous  works  we  have  presented  quantum  mechanical  treatments  of  mir¬ 
rorless  (intrinsic)  optical  bistability  (10B)  from  collections  of  spatially  dis¬ 
tributed  two-level  atoms  interacting  via  the  electromagnetic  (e.m.)  field  and 
driven  by  an  externally  applied  coherent  field. 

The  aim  of  the  present  research  project  was  to  study  the  effects  of  quan¬ 
tum  fluctuations  on  IOB  since  these  effects  have  not  been  studied  in  the 
previous  works.  The  bistability  effects  have  a  potential  significance  for  de¬ 
veloping  “optical  computers”.  The  use  of  intrinsic  bistability  as  a  switching 
mechanism  will  have  the  following  advantages: 

a)  As  the  external  feedback  such  as  mirrors  is  not  needed,  the  switching 
unit  can  be  very  small  and  therefore  much  less  complicated,  more  efficient 
and  less  expensive  than  that  used  in  a  conventional  optical  bistability. 

b)  The  switching  times  in  intrinsic  optical  bistability  will  be  short 
relative  to  the  switching  times  obtained  in  optical  bistability  with  mirrors. 
However,  the  insertion  of  noise  including  especially  the  quantum  fluctuations 
might  destroy  the  hysteresis  cycle  behavior  of  the  steady  state  response. 
Therefore,  we  have  found  it  to  be  of  utmost  importance  to  investigate  these 
statistical  fluctuations. 

Some  of  the  physical  effects  in  cooperative  resonance  fluorescence  are 
very  similar  to  those  of  the  intrinsic  bistability13-17  as  both  phenomena  are 
related  to  the  coherent  dipole-dipole  interactions.  An  interesting  problem 
in  this  field  which  has  been  studied  in  the  present  project  is  the  resonance 
fluorescence  spectra  from  two  two-level  atoms  including  the  dipole-dipole 
interactions. 
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2.  SCIENTIFIC  RESULTS 


The  detailed  scientific  results  of  the  present  project  are  given  in  the  three 
articles  added  in  the  Appendices.  Here  we  give  a  short  review  of  our  work. 

The  Langevin  equations  for  intrinsic  bistability  had  been  derived  in  pre¬ 
vious  works8-10.  The  solution  of  the  equations  was  based  on  complete  fac¬ 
torization  of  the  dipole  operators  and  then  the  Langevin  terms  drop  out  for 
the  expectation  values.  Instead  of  this  approximation  we  use  in  the  present 
work  partial  factorization,  by  which  the  expectation  values  are  made  and 
factorized  in  the  first  stage  only  relative  to  the  atoms.  The  deviation  from 
complete  factorization  is  due  to  the  Langevin  noise  terms. 

After  using  the  above  approximation  the  quantum  nature  of  the  equa¬ 
tions  is  related  to  the  “Langevin  force”  operators  /(+l(r),  which  stem  from 
the  vacuum  contribution  to  the  fluctuations,  r  is  the  retarded  time  and 
<  /‘“’Ml  > oc  S(t  -  r').  The  /(±)  operators  are  independent  of 

the  space  coordinate.  Therefore  the  treatment  of  quantum  fluctuations  in  the 
case  of  a  short  sample  (a  width  smaller  than  a  wavelength)  and  in  the  case  of 
a  long  sample  are  similar.  We  assume  that  the  /(±)  operators  are  gaussian 
operators  representing  a  gaussian  white  noise.  By  using  this  approximation 
we  transform  the  equations  for  the  inversion  of  population  and  for  the  two 
components  of  the  complex  dipole  into  three  coupled  Ito  random  differential 
equations18. 

The  quantum  statistical  treatment  of  the  intrinsic  bistability  can  be  sim¬ 
plified  for  the  case  when  the  dephasing  constant  is  much  larger  than  the 


6 


relaxation  constant  [02  >>  0i)~  In  a  pure  quantum  mechanical  model 
02  =  Pi/ 2,  but  we  may  introduce  homogeneous  broadening  by  considering 
0i  and  0i  as  empirical  cases.  In  many  experimental  cases  the  condition 
02  >>  0i  is  fulfilled.  Under  these  conditions  the  relaxation  of  the  dipole 
operators  is  very  rapid  relative  to  the  relaxation  of  the  inversion  of  popu¬ 
lation.  Therefore  for  this  case  the  expectation  value  for  the  complex  dipole 
operators  is  adiabatically  eliminated  by  using  Eq.  (17)  of  Ref.  9.  We  get 
by  this  approximation  a  one-dimensional  Ito  differential  equation  for  the  in¬ 
version  of  population.  The  methods  for  solving  such  a  random  differential 
equation  have  been  described  in  the  theories  developed  for  bistability  in  laser 
cavities19-20.  The  tunnelling  process  and  the  passage  times  between  the  two 
stable  states  in  intrinsic  bistability  were  calculated  by  using  these  methods 
for  such  cases  in  which  02  >>  0\.  The  analytical  and  numerical  calculations, 
made  by  following  this  method,  are  presented  in  Appendix  A.  Cooperative 
effects  in  the  resonance  spectra  of  two-level  atoms  are  usually  treated  by  as¬ 
suming  that  the  atoms  are  within  a  small  volume  with  dimensions  smaller 
than  a  wavelength21-22.  In  Appendix  B  we  analyze  the  resonance  fluores¬ 
cence  spectra  of  two  two-level  atoms  by  taking  into  account  the  dependence 
of  the  coherent  dipole-dipole  interactions  on  the  distance  between  the  atoms. 

In  the  present  project  we  have  treated  also  a  pure  quantum  mechanical 
system  in  which  the  relaxation  rate  0\  of  the  inversion  of  population  is  the 
spontaneous  decay  time  and  the  relaxation  rate  for  the  dipole  moment  is 
given  by  02  =  0i/2.  By  using  the  quantum  theory  of  noise  developed  by 
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Lax23  we  have  calculated  al  the  second  moments  of  the  dipole  operators  and 
studied  the  statistical  properties  of  the  reaction  field. 

It  has  been  shown  by  other  authors  that  squeezing24-26  can  be  produced 
by  optical  bistability27  but  this  possibility  has  been  studied  only  under  the 
approximation  of  a  good  cavity  limit.  For  such  cases  the  medium  relaxation 
is  much  faster  than  that  of  the  field  and  then  the  atomic  variables  can  be 
eliminated.  The  problem  of  photon  statistics  is  studied  in  Appendix  C  for  the 
extreme  case  of  a  bad  cavity  limit  where  we  preclude  any  reflecting  boundaries 
and  where  the  relaxation  of  the  field  is  much  faster  than  the  relaxation  of  the 
material.  The  statistical  properties  of  the  transmitted  radiation  are  discussed 
and  it  has  been  found  that  squeezing  is  not  obtained  in  the  present  system. 


CONCLUSIONS 

The  simple  modified  Maxwell-Bloch  equations  used  by  us  are  useful  as  a 
generic  model  for  describing  the  effects  of  quantum  fluctuations  on  intrinsic 
bistability.  The  long  passage  times  calculated  by  us,  for  a  two-level  system 
with  large  dephasing  rates  indicate  reasonably  suitable  conditions  for  obtain¬ 
ing  intrinsic  optical  bistability.  However,  the  effect  of  quantum  fluctuations 
in  terms  of  the  average  passage  times  can  be  an  important  efTect.  Quantum 
noise  effects,  as  discussed  here,  can  lead  to  significant  “glitch”  rates  which 
can  seriously  affect  stability  and  switching  properties  of  a  bistable  device. 

By  solving  our  equations  for  the  pure  quantum  mechanical  case  we  find 
that  the  reaction  field  in  intrinsic  bistability  is  bunchsd  so  that  a  quantum 
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noise  which  is  above  the  quantum  noise  limit  is  inherent  in  this  system.  We 
conclude  that  the  photon  statistics  in  the  bad  cavity  limit  is  more  chaotic 
than  the  good  cavity  limit  in  which  squeezing  was  predicted. 

The  theoretical  analysis  made  in  the  present  work  for  the  resonance  flu¬ 
orescence  spectra  of  driven  two,  two-level  atoms  can  be  used  for  evaluating 
the  spectra  in  numerical  calculations.  For  the  case  of  two,  two-level  atoms 
at  resonance  interacting  with  a  resonant  external  field,  which  are  within  a 
volume  with  a  dimension  smaller  than  a  wavelength,  we  have  derived  an  an¬ 
alytical  expression  for  the  line  shape.  Additional  side  bands  appear,  which 
are  shifted  from  the  line  center  approximately  by  twice  the  Rabi  frequency, 
but  theseifc|c bands  are  very  broad  and  very  weak. 
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l.nngcvin  nutations  for  intrinsic  optical  bistability  ( IOIJ)  in  a  system  of  two-level  atoms  interacting  by  dipole-dipole  interac¬ 
tion  arc  treated,  and  the  stochastic  clfccls  of  quantum  noise  ate  analyzed.  We  derive  a  one-dimensional  Ito  stochastic  dilictential 
equation  for  the  inversion  of  population  of  cases  in  which  the  population  relaxation  is  very  slow  relative  to  the  homogeneous 
dephasing  rate  of  the  complex  dipole.  Tor  such  eases  the  probability  distiibution,  diffusion  ptoccss  and  the  liist  aveiage  passage 
times  between  the  two  metastablc  stales  in  It  111  arc  calculated. 


Various  fot  ins  of  intrinsic  optical  bistability  ( IOB ), 
i.c.  optical  bistability  that  docs  not  depend  upon  ex¬ 
ternal  feedback  such  as  mirrors,  have  been  investi¬ 
gated  in  many  theoretical  and  experimental  works 
1 1 -I  I  ].  Thermally  induced  IOB  has  been  obtained 
experimentally  and  has  been  related  to  internal  feed¬ 
back  due  to  nonlinear  absorption  increasing  as  a 
function  of  intensity  (2,8,9).  Dagenais  and  Sharlin 
1 10)  observed  IOB  by  timing  just  below  the  absorp¬ 
tion  peak  of  the  bound  l2  cxcilon  in  CdS.  This  was 
due  to  (hctmally  induced  shift  of  the  cxcilon  reso¬ 
nance  to  low  frequency  (a  frequency  renormaliza¬ 
tion  which  is  intensity  dependent). 

Non-lhcrtnal  IOB  effects  were  observed  by  Boh- 
nert  ct  at.  (3)  and  by  Rossmann  ct  al.  (4j.  These 
experiments  were  explained  by  Schmidt  ct  al.  ( 5  J  as 
due  to  Coulomb  screening  of  carriers  caused  by  laser 
field  induced  carrier  production  resulting  in  shrink¬ 
age  of  the  baud  gap.  In  this  work  the  tuning  is  near 

1  The  icsearch  reported  herein  has  been  sponsored  in  part  by 
the  United  Stales  Army  through  the  European  Research  Of¬ 
fice  under  Contract  No.  17AM  45-87-C-0040. 


the  band  edge  in  CdS  and  is  enliiely  an  cIVcct  ol  many 
body  interaction  in  the  electron-hole  plasma. 

Currently  lull  is  routinely  observed  or  repotted  in 
the  literature.  What  has  not  been  reported  and  ap¬ 
pears  relatively  difficult  is  the  observation  of  IOB  in 
a  system  which  can  be  quantitatively  analyzed  using 
a  simple  analytical  model.  In  this  context  our  simple 
modified  Maxwell  Bloch  equations  (12-14)  aie  use¬ 
ful  as  a  generic  model.  One  should  distinguish  lie- 
tween  IOB  which  is  Ireated  in  a  steady  state  and  su- 
pci fluorescence  (SF)  (15)  which  is  a  transient 
phenomenon  depending  on  spontaneous  initiation 
anil  subsequent  coherent  pulse  buildup.  Unlike  SF, 
IOB  lias  an  externally  applied  Held  which,  when  lin¬ 
early  polarized,  causes  individual  dipole  orientation 
in  the  medium  and  the  direction  of  the  reaction  field 
of  the  atoms  is  largely  determined  by  the  incident 
driving  field,  into  the  appropriate  narrow  solid  an¬ 
gle.  Therefore,  the  dynamics  of  switching  in  IOB  is 
quite  different  from  that  of  SF  in  terms  of  initial  and 
boundary  conditions.  In  a  small  volume  (  <A‘)  SF 
cannol  happen  since  (here  is  no  preferred  diicclion. 
Whclhcr  or  nol  IOB  can  happen  for  interacting  two- 
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level  atoms  separated  by  less  than  a  resonance  wave¬ 
length,  has  been  strongly  addressed  in  the  calcula- 
tional  heuristic  treatment  of  I  Iopf  and  Uowden  [  7  J . 
In  that  work  it  has  been  shown  that  the  stochastic 
processes  in  the  two-level  atoms  with  the  dipole-di¬ 
pole  interactions  causes  bimodalily.  A  “clean  exper¬ 
iment"  which  might  be  quite  close  to  this  model, 
would  be  a  collection  of  Rydberg  atoms  (16)  in  ion 
traps  [17,18],  The  collisional  broadening  can  be 
controlled  by  laser  cooling  and  decorrelation  can  be 
controlled  by  introduction  of  thermal  IR  or  FIR  ra¬ 
diation.  Thus,  atoms  can  be  brought  closer  together 
by  laser  cooling  and  decorrelatcd  independently  by 
interacting  with  an  externally  applied  stochastic  field 
superimposed  on  Ihc  coherent  driving  field  (this 
needs  further  considerations).  As  has  been  shown  in 
the  heuristic  treatment  of  ref.  [7]  bistability  is  ob¬ 
tained  when  the  number  of  atoms  within  a  cubic 
wavelength  is  large. 

In  recent  works  1 12-14.19)  we  have  treated  the 
problem  of  IOH  in  collections  of  spatially  distributed 
two-level  atoms  interacting  via  the  electromagnetic 
Held  and  driven  by  an  externally  applied  coherent 
field.  In  the  earlier  woiks  on  IOH  of  two-level  sys¬ 
tems  the  dipole— dipole  interactions  have  been  treated 
by  assuming  a  small  volume  with  dimensions  smaller 
than  a  resonance  wavelength  [1,6]  or  by  assuming 
a  system  with  a  small  number  of  atoms  in  scmiclas- 
sical  approximation  f  7 f .  fit  (he  recent  works  [12- 
14.19)  we  have  treated  the  problem  from  the  many- 
body  standpoint  by  developing  the  Heisenberg  equa¬ 
tions  of  motion  in  the  “bad-cavity”  limit,  where  the 
relaxation  of  the  field  in  the  medium  is  much  faster 
than  the  iclaxalion  or  Ihc  material.  This  is  the  ease 
when  the  cavity  width  y  is  very  large.  Since  y " 1  cor- 
icsponds  to  the  average  photon  lifetime  in  the  cavity 
we  obtain  Ihc  extreme  "bad  cavity  limit”  by  pre¬ 
cluding  any  reflecting  boundaries  and  when  the  up¬ 
per  bound  Tor  the  photon  lifetime  is  given  by  the  pas¬ 
sage  time  through  the  medium.  If  the  cavity  width  y 
is  much  laiger  than  the  material  relaxation  rate  then 
Ihc  Held  variables  can  be  adiabalically  eliminated 
from  the  Heisenberg  equations  of  motion  with  the 
approximation  that  the  Held  is  in  steady  slate  on  the 
time  scale  of  the  material  response  time.  The  formal 
procedure  for  adiabatic  elimination,  specifically  in 
the  “bad  cavity  limit"  as  well  as  the  “good  cavity 
limit",  has  been  presented  by  Lugiato  [20],  Ihc 


“good  cavity  limit”  is  of  course  the  complement  of 
the  “bad  cavity"  approximation  and  requires  that 
the  reflectivity  boundary  condition  prevails,  i.e.  a 
high-(2  cavity.  In  this  ease  the  photon  lifetime  is  large 
and  the  cavity  width  is  small.  Thus  it  can  happen  that 
the  medium  relaxation  is  much  faster  than  that  of 
the  field  and  then  the  atomic  variables  can  be  adi¬ 
abalically  eliminated  [20,2 1  ].  An  example  would  be 
a  nonlinear  cavity  composed  of  a  Kerr  medium  be¬ 
tween  highly  reflecting  mirrors.  In  our  works  we  have 
treated  short  samples  configurations  of  two-level  at¬ 
oms  without  any  reflecting  boundary  conditions  so 
that  we  have  assumed  the  “bad  cavity  limit". 

In  our  previous  works  [  1 2-14  J  we  have  shown  that 
cubic  nonlinearilics  appearing  in  dense  two-level 
systems,  due  to  the  dipole-dipole  interactions,  can 
lead  to  IOH.  This,  of  course,  is  related  to  early  quan¬ 
tum  treatments  of  bistability  that  is  in  terms  of  fac¬ 
torization  and  cubic  nonlinearilics  1 22,23 1 .  1  lie  aim 
of  the  present  Letter  is  to  describe  the  effect  of  quan¬ 
tum  noise  on  IOB.  As  is  well  known  the  stochastic 
processes,  in  general,  change  the  phenomenon  of  bi¬ 
stability  to  bimodal  probability  distributions.  Since 
stochastic  effects  have  not  been  treated  previously 
for  IOH  our  new  analysis  of  such  effects  should  be  of 
interest. 

I  hc  equations  of  motion  obtained  for  a  short  sam¬ 
ple  of  two  level  atoms,  including  the  noise  source  op¬ 
erators,  are  given  as  [12-14] 

d(<7.(z,  t ) >„/dr  =—/?,[< a,,  T)>„  +  »] 

-I-  (,,/f,  )[/•*< a* 0(  z,  r )  >.  +  F<  „(  r,  r )  > „  ] 

-2<(7M)(z,  r )>„/'  ( r )  —  2/  '  (r)  <(7_,,(z.  t)>„ , 

(I) 

d<(7,„(z,  r)>.-,/dr=  -  (/<£/2/t)  <o.(z,  r)>„ 

+  i[/j-c<o,(z,  r) >„ j <rTf(t(z,  r)>, 

-/?2<at„(z,r)>.+/-(r)  <eif(r,  t)>, .  (2) 

Here  (a,n)„  is  the  expectation  value  of  the  slowly 
varying  complex  atomic  polarization  per  unit  vol¬ 
ume,  and  <or>„  is  the  expectation  value  of  the  atomic 
inversion  per  unit  volume.  The  average  <  )„  are  with 
respect  to  the  atomic  slates  only.  The  parameter  /t  is 
the  modulus  or  the  matrix  element  of  the  transition 
dipole  moment  or  an  atom,  and  n  is  the  number  of 
atoms  per  unit  volume.  d=io—io„  is  the  detuning 
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where  to  is  the  frequency  of  the  applied  field  E  and 
co„  is  the  atomic  resonance  frequency.  The  rales  p[ 
and  /?2  are  the  inverses  of  the  relaxation  times  7',  and 
T2  for  the  inversion  of  population  and  the  dipole 
moment,  respectively,  and  d=m-ro0  is  the  devia¬ 
tion  of  the  applied  field  frequency  to  from  the  atomic 
resonance  frequency  <o„.  t  is  the  frequency  renor¬ 
malization  constant  derived  in  our  previous  work 

t  =  Tni)flc3/4to*i,  (3) 

where  p  is  the  spontaneous  decay  constant 

P=4ft2<Oo/3fH'* .  (4) 

The  frequency  renormalization  stems  from  coherent 
dipole-dipole  interactions  between  atoms  that  are 
within  a  volume  of  a  cubic  wavelength.  The  electro¬ 
magnetic  lic-lil  appeals  in  cqs,  (  I  ),  (2)  in  two  pails, 
i.e.  the  free  field  /  and  the  macroscopic  field  E.  1  lie 
latter  is  assumed  to  be  a  constant  classical  driving 
field  (or  equivalently  that  il  is  in  coherent  stale).  The 
free  field  operators  of  eqs.  ( I  ),  ( 2 )  obey  in  the  plane 
wave  limit,  the  equations 

<f  ,(r)/  (r')>=/W(T— r  )  , 

(f  MO./'  ‘  (r)>  =  </-(r)/-(r')> 

=  </-(r)/'(r')>=0,  (5) 

where  the  averages  <  >  arc  expectation  values,  and 
r  is  the  retarded  time  variable  t  =  /-z/c.  These  op¬ 
erators  constitute  fluctuating  noise  terms  which  stem 
from  the  quantum  field  coupling  to  the  atoms. 

In  our  previous  works  we  have  treated  the  prob¬ 
lem  of  factorization  of  products  of  atomic  operators 
that  arise  from  adiabatic  elimination  of  the  reaction 
field  contributions.  Here  the  multiplications  with  the 
slowly  varying  external  field  E  are  automatically  fac¬ 
torized  under  the  present  assumption  of  a  driving 
classical  field.  The  only  operator  product  terms  left 
in  cq.  ( I )  are  the  multiplicative  noise  terms,  anil  this 
of  course  factorize  with  respect  to  the  atomic  initial 
states  (in  consistence  with  the  markovian  approxi¬ 
mation  of  eq.  (5) ).  In  the  present  work  we  assume 
also  adiabatic  elimination  conditions  for 
which  means  overwhelming  large  dephasing  of 
<<7»o>  (Pi2** Pi )  Thus,  under  these  conditions,  the 
coherence  of  the  off-diagonal  dipole  terms  is 
quenched  and  the  diagonal  and  off-diagonal  terms 
are  dccorrclatcd.  Thus,  under  the  adiabatic  approx¬ 


imation  for  strong  dephasing,  deem  relation  or  cj,„ 
and  oz  is  an  excellent  approximation,  and  indeed, 
must  follow  for  consistency.  One  should  be  aware  of 
the  fact  that  the  effect  of  lOll  treated  here  is  not  re¬ 
lated  at  all  to  the  cooperative  effects  discussed  for 
SI7.  It  is  related  to  the  renormalization  of  the  fre¬ 
quency  due  to  the  high  density  of  the  two  level  atoms 
and  is  similar  to  the  renormalization  of  frequency 
for  IOII  in  semiconductors  as  explained  in  the  in¬ 
troduction  to  this  Letter. 

By  adiabatic  elimination  of  < rr , ,,  >  and  (n  „>  in 
cqs.  (1),  (2)  we  obtain  the  single  Ito  equation 
[24,25]  for  the  inversion 


di 7 


/?.("-  +  I) 


-(■ 

(  -<Uu  a. 

\i(,l  — 1(7. 


II 


-  d  /  *  +  - 

)  l(d  —  to 


|W|,  IV/L  \  . 

-T.-)7/i.r' 
fi?i-  d/7  )  . 

<7)1/1-  / 


(6) 

where  (oK  =  /(/.'/ ti  is  the  Rabi  licquency  of  the  ap¬ 
plied  external  field.  Here,  we  have  defined  the  di¬ 
mensionless  stochastic  variables  d /  1  and  d/7  by 
the  equations 


d£'(r) 

dr 


=/  '  (r). 


d/Jj) 

dr 


(r)  , 


(7) 


and  have  dropped  the  <  >.,  notation  for  convenience 
and  scaled  the  atomic  variables  to  functions  corre¬ 
sponding  to  densities  per  atom.  The  stochastic  pmp 
ci  ties  of  the  variables  d/-‘ '  and  d /  follow  from  cqs. 
(5) and  (7), 

<dF  *  ( r )  d /•'  (r  )>=/?  cl r , 

<dF'(t)d/-  '  (r  )>  =  <d/-  ( r)  <1/  (r')> 

=  <d/'  (r)dF'(t'  )>=().  (8) 

T  he  mode!  can  be  easily  generalized  to  other  spe¬ 
cific  radiation  rose;  voir  conditions,  such  as  a  ther¬ 
mal  reservoir  at  temperature  /VO.  Using  conditions 
(8)  we  transform  cq.  (6)  into  a  Fokker-I’lanck 
equation  [24,25]  for  the  inversion  of  population  o: 
as  a  function  of  the  parameters  |ruK|J,  J,  f,  //,  and 
Pi  ( where  /?2 »  ft,): 
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jj/lffr) 

Or 


J)_ 

i)(7. 


[( 


-fl,{n2+l) 


l"'i.  \'n  ft-  ,  .1 

1  0*  /  IwrI-o:/?!  .  7 

2ito?  \P\  +  (A-«t:YiJ(  :)J-  (9) 


We  represent  the  drift  term  /l(o;)  and  the  diffusion 
cocllicicnt  B ( a. )  as 


.  I  ( rr. )  =  -  /?,  ( cr-  +  1  )  - 


\«>k\2OzPi 
ft\  +  (A-tn.)2  ’ 


/f  ( rr_ )  = 


_  |wK  I  ’oj/?, 
/?!+" "(d-f(7r)5  ■ 


(10) 

(ii) 


The  steady  slate  solution  of  (9)  is  given  by 


/ M,Vr)  =  C'j;;— rexp  2 

B{n.) 


j  -/?(»-)/ 


r 

Bin  ) 


exp(-l/(ff.)|  , 


(12) 


wlicte  (is  a  normalization  constant  and  the  poten¬ 
tial  V(o.)  is  defined  as 


t'in:)=- 2 


.tin. 

Bin,)' 


(13) 


i’in.)  plays  a  role  similar  to  a  chemical  potential 
since  the  main  dependence  on  a.  is  through  this 
function  while  B{ar)  is  a  slowly  varying  function  of 

17.. 

Hy  a  straightforwaid  integration  we  gel 


f'('7,)=ln(  \n.  | ) 


"2/7, 

j; + 


2(ft]  +  J?)-4JA 
|wR|J  ) 


+  (1+0.) 


2v  -A/k 
I  Wr  1 3 


+  (o:-l) 


7 

e 

k'7nlJ 


Mlti'A 
i  (nH  | J 


(14) 


I o  illustrate  the  quantum  noise  effects  in  intrinsic 
bistability  we  show  in  fig.  I  the  value  of  <<r,>  as  a 
function  of  |w„|J,  where  all  rates  are  normalized  to 
ft,.  T  his  figure  has  been  calculated  according  to  the 
steady  state  cubic  equation  of  stale  obtained  in  our 


t  ig.  I.  (a,)  as  a  function  of  pu„|!  Tor  normalired  rotameters 
/?,  =  I.  /?,=  10,  4=  -  10  and  r  =  100,  wltctc  all  quantities  arc  in 
units  of  ft,,  the  solid  and  dotted  curves  represent  respectively 
stable  and  unstable  steady  stales.  The  points  a.  b  and  c  represent 
the  lliice  solutions  of  the  cubic  steady  stale  equation  where  the 
subscripts  7=1,  2  and  3  refer,  respectively,  to  external  fields 
| ii)K  | J  =  1 00,  120  and  1 50  in  units  of/1,. 


pteviotts  work  (ref.  [12],  cq.  (18)).  and  this,  of 
course,  is  in  total  agreement  with  the  expression  for 
the  drift  term  of  eq.  (10).  T  he  dephasing  rate,  ft,, 
has  been  chosen  to  be  10  ft,,  so  that  the  approxi¬ 
mation  /?2 »  ft,  is  justified.  To  obtain  a  condition  of 
intrinsic  bistability  the  parameter  t  must  satisfy 
[6,14]  f//??>6.  For  illustration  we  have  chosen 
t=  100  and  A-  -  10,  in  units  of  ft,.  T  his  gives  the 
ratio  c//?2=  10.  The  solid  and  dotted  curves  in  fig.  I 
represent,  respectively,  stable  and  unstable  steady 
states.  T  he  points  a,  b  and  c  in  this  figure  represent, 
respectively,  the  three  steady  state  solutions  Tor  a. 
for  a  prescribed  Held  value  |«7„|\  where /=  I,  2,  3 
refer  to  | o>l(  | 2  =  100,  120,  150  respectively.  For  the 
same  conditions  and  prescribed  fields  the  corre¬ 
sponding  probability  functions,/""  (o.),  calculated 
from  eq.  (12),  arc  depicted  as  functions  of  n.  in  fig. 
2.  The  maxima  of  the  probability  density  functions, 
in  each  ease,  correspond  to  the  two  stable  solutions, 
a  and  c,  or  the  cubic  equation  whereas  the  minima, 
b,  correspond  to  the  associated  unstable  solution  of 
the  equation. 

Hie  time  dependent  Fokkcr-I’lanck  equation,  cq. 
(9),  can  be  described  by  [20] 

T)/i)t  =  A/(oI,  t)  , 


(15) 
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Fig  2.  The  probability  density  /(a.)  for  various  external  fields 
as  a  function  of  n:  for  stationary  states  with  normalized  paramc- 
leis/J,  =  \,Pj  =  10,  ./  =  -  10, « =  100. 1  lie  three  curves  arc  for  Fields 
with  |ojk|!=I00,  120  and  ISO.  The  maxima  a  and  c  correspond 
to  the  stable  solutions  of  the  cubic  equations  and  the  minima  h 
to  the  corresponding  unstable  solutions,  where j=  I,  2,  3  refer  to 
|(uk|3=  100,  1 20 and  1 50,  respectively.  All  paramcleisarc  in  units 

or/),. 

where,  under  conditions  similar  to  those  which  per¬ 
tain  to  llg.  2,  the  linear  operator  exhibits  a  nearly  de¬ 
generate  pair  of  lowest  lying  eigenvalues,  A,,  =  0  and 
A ,  A  0,  and  a  large  gap  between  A ,  and  the  remaining 
eigenvalues.  The  lime  evolution  off  (a.,  r)  occurs  in 
two  separate  steps.  In  the  first  one  the  relaxation 
process  is  local  in  which  the  two  peaks  of  the  double 
peaked  distribution  described  in  fig.  2,  evolve  in¬ 
dependently  of  each  other  (20).  This  stage  is  very 
rapid  corresponding  to  the  high  eigenvalues  for  which 
A(»A|,  (j>  I ).  The  shape  of  the  two  peaks,  of  a  given 
distribution  at  the  end  of  the  first  stage  is  identical 
to  that  represented  by /“"(a.),  but  the  ratio  of  the 
maxima  of  the  peaks  may  be  different,  depending  on 
the  initial  conditions.  Referring  to  fig.  I ,  we  may  start 


in  the  upper  (highly  excited )  branch  of  a.  and  change 
|  toK  | 2  in  the  neighborhood  of  a  to  approach  the  point 
a,  or  we  may  start  on  the  lower  branch  and  arrive  at 
the  point  c.  After  the  fast  relaxation  of  the  first  step 
the  probability  density/ ( rr- )  is  given  by  the  line  pro¬ 
file  about  the  point  a,  whereas  for  the  second  case  it 
is  given  by  the  line  profile  about  point  c,  as  illus¬ 
trated  in  fig.  2. 

In  the  long  time  limit  (the  second  stage)  the  func¬ 
tion  J(a2)  can  be  approximated  by 

f(az,  r)=/<s,,((7r)  +  cxp( -A,  r)/,(<7-)  ,  (16) 

where  /,  ( crr )  is  the  eigenfunction  of  the  Fokker- 
l’lanck  operator  corresponding  to  the  eigenvalue  At. 
The  probability  to  find  the  system  in  the  slates  a  and 
c  at  lime  r  arc  given,  respectively,  by 
0 

"’»(0=  J  /(o_.)  da. , 

IV,  (  t  )  =  I  -  tV,  ( t )  .  (17) 

The  rate  by  which  iv.,(r)  changes  is  given  by 

dMv(r)/dr=-(l/T,)iv.  +  (l/r,.)(l-»r.)  ,  (18) 

where  ra  and  rt.  arc  the  “escape  times”  out  of  the  a 
and  c  wells  of  the  potential.  On  the  other  hand  we 
have 

dit'„/dr=  —  Al|tvil(r)  — ii\,(co)  )  .  (19) 

By  comparing  cqs.  (18)  and  (19)  we  gel 

(20) 

By  using  methods  described  in  refs.  [  20,21  |  and  ap¬ 
proximations  that  are  valid  within  the  switching  re¬ 
gime,  the  “escape  times”  are  given  by  |20| 


0 


ta=/W  [  da.  [A/(a.)  j 

|-'cxp[-(/(o.)| 

n;{b) 

r 

tc  =M  J  da,  [11(0,)] 
-  1 

- 1  cxp[  -  (/(a.  ) ) 

where  M  is  a  function  or  the  potential  barrier  be¬ 
tween  the  two  potential  wells. 
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iU  =  2  |  dfx.  exp |  V{a.)  ]  .  (22) 

( 

By  calculating  the  escape  times  according  to  the 
above  equations  for  the  Ihrcc  sets  of  conditions  il¬ 
lustrated  in  ligs,  I,  2,  we  obtain  r(u2) % 3. Ox  I01 
P\  r(rj)«l.65x  10’  ft1-,  t(a,)*92.6  ft1, 
r(r, ) « 3.08 x  10'  r(<i,)»5.6x  10'  ft1, 

r ( ft )  =  36.7/?j  l  or  the  probability  density  corre¬ 
sponding  to  |o7„  |  ‘  =  120,  labelled  j=  2  in  fig.  2,  we 
obtain  a  bimodal  distribution  and  long  passage  times. 

I  bis  condition  corresponds  very  nearly  to  a  gener¬ 
alized  Maxwell's  construction  and  the  calculated 
passage  times  arc  very  close  to  the  maximum  aver¬ 
age  passage  time  for  the  parameters  chosen.  In  re- 
gaid  to  the  other  two  eases  illustrated  in  figs.  I,  2  for 
/=  1 .  3  corresponding  to  |  r»„  | 2  =  I  (X),  150  respec¬ 
tively,  the  probability  densities  aic  nearly  monom- 
odal  with  drastically  disproportionate  escape  times 
corresponding  to  the  two  stable  steady  slate  regions. 

I  lie  long  passage  times  associated  with  the  regime 
of  dearly  distinct  bimodality  indicate  reasonably 
suitable  conditions  for  observation  of  IOB  in  prin¬ 
ciple.  but  the  effect  of  quantum  fluctuations  is  a 
measurable  quantity  in  terms  of  the  average  first 
passage  times  anil  ran  be  an  important  effect,  es¬ 
pecially  in  the  appioach  to  ncaily  mouomoda!  be¬ 
havior  near  the  steady  state  switching  thresholds. 
Quantum  noise  effects,  as  discussed  here,  can  lead  to 
significant  “glitch”  rates  which  can  seriously  affect 
stability  and  switching  properties  of  a  bistable  device. 
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Resonance  Fluorescence  Spectra  of  Two  Driven  Two- 

Level  Atoms 

YACOB  BEN-ARYEJi  and  CHARLES  M.  BOWDEN 


Abstract— We  develop  a  theoretical  analysis  for  the  spectrum  oT  the 
light  scattered  by  two,  two-level  atoms  driven  near  resonance  by  an 
external  coherent  electromagnetic  licit!.  The  atoms  are  assumed  to  re¬ 
lax  to  equilibrium  with  the  driving  field  due  to  radiation  damping.  The 
explicit  dependence  of  (he  coherent  dipole-dipole  interactions  on  the 
distance  between  the  two  atoms  Is  taken  into  account  in  the  theoretical 
treatment  of  the  spectra.  For  the  case  of  two,  two-level  atoms  with  a 
fixed  distance  smaller  than  a  wavelength,  an  analytical  result  for  the 
resonance  fluorescence  line  shape  Is  obtained,  and  the  result  Is  com¬ 
pared  to  previous  works. 

1.  Introduction 

NE  of  (lie  most  interesting  subjects  in  quantum  optics 
which  has  been  studied  extensively  in  the  litcraluic 
is  the  resonance  fluorescence  spectra  from  a  two-level 
atom  interacting  with  a  nearly  resonant  coherent  external 
field  |l]-[5|.  Interesting  quantum  effects,  such  as  photon 
antibunching,  have  been  observed  |6|  in  connection  with 
the  resonance  fluorescence  phenomena.  The  structure  of 
the  spectrum  of  resonance  fluorescence  of  a  cooperative 
system  of  atoms  has  been  investigated  in  many  woiks  |7|  - 
1 1 3 1 .  Scnitzky  |7|  was  the  first  to  predict  that  the  spectrum 
of  a  cooperative  system  would  consist  not  only  of  the  usual 
peaks,  but  also  of  additional  side  bands  at  the  harmonics 
of  the  Rabi  frequency.  This  problem  has  been  investi¬ 
gated  further  in  the  various  works  |8]-|13|,  and  in  partic¬ 
ular,  the  dynamics  of  the  two,  two-level  atoms  lias  been 
investigated. 

The  purpose  of  the  present  work  is  to  study  resonance 
fluorescence  spectra  of  two,  two-level  atoms  by  the  use  of 
new  methods  and  compare  the  results  to  other  works. 

The  spectral  distribution  of  the  resonance  fluorescence 
scattered  light  of  one  two-level  atom  is  well  known  |l|. 
It  has  been  well  established  experimentally  |I4|  that  for 
low  incident  field  strengths,  the  scattered  spcctium  con¬ 
sists  of  a  single  broadened  line,  which  above  certain 
threshold  strength,  it  exhibits  a  pair  oT  side  bands  in  ad¬ 
dition  to  the  central  component.  The  ratio  of  the  central 
line  to  the  side  band  peak  heights  is  3:1,  while  the  line 
width  ratio  is  1 :  1 .5.  This  kind  of  spectrum  is  in  complete 
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agreement  with  that  of  Mollow  [I],  Oliver  el  al.  1 1 5 J , 
Carmichael  and  Walls  [4),  Kimble  ct  al.  |6J,  and  many 
others.  This  result  is  in  contrast  with  the  one-photon  ap¬ 
proximation  1 16],  which  gives  a  ratio  2  :  1  rather  than  3  :  I 
in  the  peak  heights.  A  more  sophisticated  one-photon 
model  for  resonance  fluorescence  spectra  of  a  two-level 
atom  has  been  developed  by  Baklanov  [17]  and  by  Sten- 
liolm  1 18],  This  new  one-photon  model  gives  the  correct 
resonance  fluorescence  spectra  and  is  in  complete  agree¬ 
ment  with  other  works  11],  |4|,  [6],  ]I5].  It  is  the  purpose 
of  the  present  work  to  apply  this  approach  1 1 7] ,  [18]  for 
developing  the  theoretical  analysis  of  the  resonance  fluo- 
icsccncc  spectra  of  two,  two-level  atoms  and  to  compare 
the  results  to  those  of  other  authors  [8]~|I3]. 

Cooperative  atomic  effects  in  the  resonance  fluores¬ 
cence  spectra  arc  usually  treated  by  assuming  that  the  at¬ 
oms  arc  within  a  small  volume  with  dimensions  smaller 
than  a  wavelength.  The  nature  of  the  coherent  dipole-di- 
polc  interactions  between  two  identical  two-level  atoms 
at  a  fixed  arbitrary  distance  r  apart  lias  been  studied  by 
Milonni  and  Knight  |I9|.  The  coherent  dipole-dipole  in¬ 
teractions  with  their  explicit  spatial  dependence  may  ex¬ 
plain  various  cooperative  effects,  such  as  supcrradiance 
[20]  and  intrinsic  bistability  [21].  In  the  present  work,  we 
analyze  the  resonance  fluorescence  spectra  by  taking  into 
account  the  dependence  of  the  coherent  dipole-dipole  in¬ 
teractions  on  the  distance  r  between  the  atoms. 

In  Section  II,  we  present  the  general  theory  for  the  dy¬ 
namics  of  driven,  spatially  separated  two,  two-level  at¬ 
oms.  We  follow  here  the  approach  given  in  [19],  by  which 
the  dipole-dipole  interactions  lead  to  a  damping  mecha¬ 
nism,  but  in  addition,  we  include  here  the  effect  oT  the 
coherent  interaction  with  the  externally  applied  field. 

In  Section  III,  we  develop  the  density  matrix  equations 
of  motion  by  generalizing  the  calculations  made  in  |I7| 
and  [18]  for  one  atom  to  two,  two-level  atoms.  The  anal¬ 
ysis  for  the  resonance  fluorescence  spectra  is  made  in  Sec¬ 
tion  IV.  Although  for  the  general  cases  it  was  not  possible 
for  us  to  obtain  final  analytical  results,  the  equations  de¬ 
veloped  by  us  can  be  used  further  in  a  numerical  analysis. 
Also,  wc  have  emphasized  the  analysis  of  specific  cases 
for  which  final  physical  conclusions  have  been  made. 

II.  fjENDRAL  Theory 

Our  system  is  composed  of  two  identical  isolated  two- 
level  atoms  fixed  in  position  with  a  distance  r  apart,  driven 
by  a  strong  prescribed  incident  coherent  electromagnetic 
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field,  which  is  nearly  resonant  with  the  two-level  atomic 
frequency.  It  is  assumed  that  the  two  atoms  come  into 
equilibrium  with  the  driving  field  through  the  effect  of  ra¬ 
diation  damping. 

The  Hamiltonian,  which  describes  (he  system,  is  given 
in  the  rotating  wave  and  in  the  electric  dipole  approxi¬ 
mation  by  //  -  //„  +  //'  +  H": 

H0  =  E  <t'°  +  ft  £  uqaqaq 

/  l  -  I  q 

2 

//'  =  -ih  E  2  +  lie 

'=1  q 

-7  2 

II"  =  E  uRo'iV *“•  ?,)  +  HC  (2. 1 ) 
where  1 1C  denotes  the  Hcrmitial  conjugate  and 


^  \~T, 7 

Hn  includes  the  free  atoms  and  free  field  Hamiltonian.  We 
use  here  the  usual  operators  related  to  the  SU(2)  algebra 
fora  single  atom:  ctJ'1  represents  (lie  population  inversion 
operator  and  Oj1  are  the  raising  and  lowering  operators 
for  atom  i  with  coordinates  r,.  wR  is  the  Rabi  rate  asso¬ 
ciated  with  the  applied  coherent  field.  w0  and  k„  are  the 
carrier  frequency  and  wave  vector  of  the  applied  field, 
respectively.  V  is  the  quantization  volume  for  the  field, 
oiq  is  the  frequency  and  eq  is  a  unit  vector  in  the  direction 
of  polarization  for  the  mode  q,  jt(  is  the  dipole  moment 
vector  of  atom  f,  and  o>  is  the  atomic  frequency  of  the  two- 
level  system. 

We  limit  the  discussion  in  the  present  paper  to  the  case 
of  an  external  field  propagating  perpendicular  to  the  in¬ 
teratomic  axis.  By  choosing  the  center  of  coordinates  at 
(he  center  of  the  interatomic  axis,  one  finds  immediately 
that  the  spatial  factors  of  //"  vanish.  Therefore,  for  this 
case,  the  interaction  with  the  radiation  is  completely  sym¬ 
metric  relative  to  the  two  atoms  and  the  analysis  of  the 
resonance  fluorescence  spectra  is  relatively  simpler.  The 
treatment  of  the  general  case  where  the  radiation  is  prop¬ 
agating  in  any  direction  can  be  done  in  a  straightforward 
way  by  following  the  present  methods,  although  the 
expressions  will  be  more  complicated  (the  interaction  II" 
inserts  spatial  factors  in  the  matrix  elements  which  are 
different  for  the  two  atoms),  for  cases  for  which  k„  ■  r 
«  I,  the  spatial  factors  of  //"  can  be  neglected  and  the 
present  analysis  remains  valid  for  any  direction  of  the  ex¬ 
ternal  applied  field. 

A  complete  set  of  the  atomic  states  can  be  given  as 

k>  -!  +  >,!  +  >,; 

=  ^[|  +  ),|->2  +  i->.i  +  >2!; 

l*3>  =^II  +  >,|->J  -  I  -  >,l  +  >J; 

I*)  -  I ->■!->,  (2,3) 


where  |  T  >  and  |  -  >  refer  to  an  atom  in  the  upper  and 
lower  state  of  the  two-level,  respectively,  and  the  indexes 
I  and  2  refer  to  the  atom  number. 

In  the  Hamiltonian  of  (2.1),  we  distinguish  between  the 
coherent  interaction  //"  with  the  external  applied  field  and 
the  incoherent  interaction  H'  with  the  continuum  of  the 
electromagnetic  modes.  The  incoherent  interaction  de¬ 
scribes  the  escape  of  excitation  energy  into  the  radiation 
field  and  is  noted  as  a  dissipative  mechanism.  This  dis¬ 
sipation  appears  in  a  way  that  is  analogous  to  a  relaxation 
to  a  heat  bath  where  the  heat  bath  is  understood  as  the 
quantum  fluctuations  of  the  photon  vacuum. 

The  influence  of  the  incoherent  interaction  //'  on  the 
equations  of  motion  can  be  calculated  separately  from  that 
of  the  coherent  interaction  in  the  following  way.  By  as¬ 
suming  that  the  atomic  system  is  in  initial  state  |  tj>,„  )  (m 
=  I  •  •  -4),  its  decay  into  a  lower  state  is  evaluated  by 
developing  the  equations  of  motion  due  to  II'  and  solving 
them  by  the  use  of  the  Laplace  transform  within  the  pole 
approximation  (see  explicit  calculations  in  1 19)). 

We  find  that  the  damping  constants  for  the  states  |4>2  > 
and  |  (pi  >  are  given  by  /3+  and  _ ,  respectively,  where 

fit  =  0(  I  ±  *(r))  =  0  ±  B;  B  =  fig(r).  (2.4) 

Here  0  =  2pJw5/3/»c1  is  one  half  of  the  Einstein  A  coef¬ 
ficient  and  g(r)  is  a  spatial  factor  which  depends  on  the 
distance  between  the  atoms.  g(r)  is  described  explicitly 
in  |I9|  or,  in  another  form,  B  =  0g(r)  is  described  in 
(21).  g(r)  tends  to  I  forr  -*  0  and  tends  to  0  for  r  -*  °°. 
The  imaginary  part  of  0g(r)  =  B  corresponds  to  fre¬ 
quency  shifts  and  the  absolute  values  of  the  amplitudes  of 
|  <t>i  >  and  |  >  decay  at  a  rate 

Re  /?£=/?(  I  ±  Re  g(r))  =  0  ±  Re  B.  (2.5) 

In  the  same  way,  we  find  that  the  amplitude  for  the  state 
|  '/>(  >  has  a  damping  constant  20  where  the  decay  constant 
from  the  state  |tf>i>  to  the  slate  |  02  )  is  0,  =  0(\  + 
g(r))  =  0  -I  B,  while  that  front  |<f>t>  to  |  <#>3  >  is  0 -  = 
0(\  -  g(r))  =  0  —  B.  Taking  into  account  the  imagi¬ 
nary  pail  of  g(r)  (including  frequency  shifts),  the  effect 
of  //’  on  the  amplitudes  of  the  states  |  0„,  >  can  be  de¬ 
scribed  by  the  following  diagram. 

|<M 


In  the  present  work,  we  calculate  the  resonance  fluo¬ 
rescence  spectra  in  a  strong  field  by  including  in  our  re¬ 
duced  system  two  sets  of  stales:  the  states  (m  = 
I  •  4)  which  do  not  include  scattered  photons,  and  the 

states  |  <pm  q  >  which  include  one  scattered  photon  in  the 
mode  q.  We  introduce  the  atomic  reduced  density  matrix 
elements 

Pma  ~  (0«|p«.UB>;  Pin. mi  =  <  <t>m  | Pol  I  't’n.q  )  (2-6) 


1.178 


where  we  have  traced  out  the  slates  of  the  continuum  held. 
We  calculate  the  equations  of  motion  for  these  density 
matrix  elements  by  using  the  Hamiltonian  //„  +  //"  of 
(2.1).  The  effect  of  //'  on  the  equations  of  motion  has 
been  reduced  to  the  use  of  the  damping  constants  evalu¬ 
ated  previously  for  the  amplitudes  of  the  states  |  <f>,„  >  (or 
Under  steady-state  conditions,  the  calculation  of 
the  density  matrix  elements  p,„  „  and  p,„„n  is  reduced  to 
the  solution  of  algebraic  equations.  The  procedure  is  com¬ 
pletely  analogous  to  the  method  used  in  (I7|,  (18).  How¬ 
ever,  the  calculations  for  our  system  of  two,  two-level  at¬ 
oms  become  much  more  complicated  than  those  made  for 
one  two-level  atom.  The  complication  is  mainly  due  to 
(he  fact  that  each  of  the  sets  p,„  „  and  p include  16 
matrix  elements  (m  =  I  •  •  •  4,  n  =  I  •  •  ■  4)  and  the  so¬ 
lution  of  16  simultaneous  algebraic  equations  becomes 
quite  tedious.  Although  in  the  present  paper  we  develop 
the  theory  for  general  eases,  we  emphasize  finally  special 
eases  for  which  the  calculations  and  the  physical  solutions 
become  quite  simple. 

The  theory  and  the  calculations  which  arc  presented 
here  can  be  used  further  for  evaluating  the  resonance  flu¬ 
orescence  spectra  of  two-level  atoms  by  numerical  cal¬ 
culations  solving  the  algebraic  equations.  In  the  following 
section,  we  present  the  calculations  for  the  reduced  den¬ 
sity  matrix  elements  p,„  „  and  p,„  ,,n.  These  matrix  elements 
arc  used  for  evaluating  the  resonance  fluorescence  spectra 
in  Section  IV. 
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p22,  and  use  instead  the  normalization  condition 

Pi.  +  P22  +  P33  +  P44  =  I-  (3.5) 

By  assuming  steady-slate  conditions,  i.e.,  p„„,  =  0,  we 
get  16  equations  for  the  density  matrix  elements  p„,„  (in, 
n  =  I,  •  •  •  ,  4).  All  the  density  matrix  elements,  which 
are  coupled  to  the  antisymmetric  slate  |  <#>3  > ,  can  be  elim¬ 
inated  from  this  equation  in  the  following  way. 

From  the  equation 

Pi,  =  0  =  2(0  -  Re  B)Pn  -  2(0  -  Re  B)p„  (3.6) 
we  get  the  result 

Pi: i  —  Pi  i-  (3-7) 

The  set  of  equations 

Pn  =  0  =  -iApij  -  inpu  ~  (30  -  D*)p,\ 

Pu  =  0  =  —inpty  ~  iahx  ~  2(0  -  Im  B)pn 

P.14  =  0  =  -TAp34  +  1C*PJ2  -  (0  -  fl)P34  (3.8) 

(including  also  their  complex  conjugates)  are  homoge¬ 

neous  algebraic  equations  whose  only  solutions  are 

Pti  =  P23  =  Pl4  =  0-  (3.9) 

The  equations  for  the  density  matrix  elements,  which 
arc  not  coupled  to  the  state  |  ) ,  become 

Pi.  =  0  =  icr(pi2  -  P2I)  -  40pn; 


111.  This  Disnsitv  Matrix  liyttAiioN  or  Moiion 

We  introduce  the  slowly  varying  nondiagonal  density 
matrix  elements  by  using  the  definitions 

Pn  =  P-2P'W  =  Pn:  Pu  =  Pn^"'  =  (/»«)* 

Pi.i  =  Pti^"  =  (P,,)*;  hi  =  P.i4P'w  -  (hi)* 

P.4  =  Pi  4CJ,W  =  (P4,)*l  P23  =  P23  =  (hi)*  (3.1  ) 

where  the  deviation  from  resonance  is  given  by 

A  =  uj  —  too .  (3.2) 

We  choose  the  phase  of  (oR  so  that  it  is  an  imaginary  num¬ 
ber  and  define 

ww  =  (a  V2  (3.3) 

where  a  is  a  real  number.  This  choice  of  phase  corre¬ 
sponds  to  the  Hamiltonian  of  interaction  and  the  interac¬ 
tion  constant  a  defined  in  [  1 8 1  for  treating  the  resonance 
fluorescence  spectra  for  one  two-level  atom. 

In  addition  to  the  damping  of  the  density  matrix  ele¬ 
ments  p„m,  we  include  in  these  equations  the  increase  of 
the  population  of  the  levels  |<£2>,  and  |<f>4),  rep¬ 

resented  by  p22,  pji,  and  p44,  due  to  the  decay  from  the 
upper  states.  The  equations  for  (he  diagonal  matrix  ele¬ 
ments  arc  dependent  via  the  relation 

Pit  +  P22  +  P.13  +  P44  =  0  (3.4) 

so  that  we  should  eliminate  one  of  the  equations,  e  g.,  for 


2p.i  T  Pn  +  P44  =  I 

P44  =  0  =  icy ( p42  -  P24)  •  2(0  -I  Re  B)pn 
T  2(0  -  Re  B)pn 

P 14  =  0  =  -2iApl4  +  i«(p,2  -  p24)  -  2 0p,4 


P12  =  -rApl2  +  i«[pn  -  pn  +  Pm]  -  (30  +  B*)pn 

Pu  =  0  =  ~iAp24  +  icy[p22  -  P44  +  Pm] 

-(0  +  fl)P24  (3.10) 

where  the  equations  for  pl2,  p24,  and  p )4  are  complex, 
while  other  c(|uatious  arc  real.  These  equations  may  be 
easily  calculated  by  numerical  analysis. 

For  a  very  strong  external  field  at  exact  resonance,  the 
approximate  solutions  of  (3.10)  arc  given  by 

Pn  —  P22  =  P33  ~  P44  -  7 1  Pm  =  0 
4 


P21 


P42 


2i0Ph 

Of 


(for  A  =  0; 


cy  »  0,  a  »  |0|). 


(3.11) 


Equations  (3.11)  represent  the  fact  that  when  the  pump¬ 
ing  by  the  driving  field  at  resonance  (A  =  0)  is  very 
strong,  the  population  in  the  four  stales  becomes  approx¬ 
imately  equal.  In  the  present  paper,  we  discuss  analyti¬ 
cally  the  problem  only  for  this  limiting  case.  For  other 
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cases,  we  represent  the  solutions  by  deriving  sets  of  al¬ 
gebraic  equations  which  can  be  solved  numerically. 

We  introduce  slowly  varying  density  matrix  elements 
p,„„,  using  the  following  definitions: 

Pl.iq  ~  Pi  lqe  1"*"<  Pi.  Iq  Pl.iq'  Pi,  Iq  =  Pi.  Iq?  3l“'“ 


Pl.iq  ~  Pl.iq? 


Pl.iq  ~  Pi.  iq'  Pi. iq  ~  Pl.t, I1 


Pi. Iq  —  Pl.iqP'^';  Pl.iq  —  Pl.iq?  Pl.iq  ~  Pl.iq? 

Pl.iq  ~  Pl.iq''  Pl.iq  =  Pl.iq?  Pl.iq  =  Pi ,  iq  ?  ~  2' 

(3.12) 

where  w0  is  the  frequency  of  the  external  field,  l;or  these 
density  matrix  elements,  we  get  under  steady  conditions, 
three  independent  sets  of  algebraic  equations. 

For  p,„  ,,,  where  m  #  3,  we  get  nine  algebraic  equa¬ 
tions: 

Pl.iq  =  0  =  ->lqp II  f  (  -i(&  -  6)  -  3/3  -  n'  \f>t.lq 

1  ,nlPl.lq  ~  Pl.iq  3  Pi  .iq  | 

Pi,\q  =  o  =  [/'( 2A  h  5)  -  2  /3  ]/?.,,,  I  ini  Pl.iq  ~  Pi.  ij 
Pl.iq  =  0  =  [/(A  6)  -  3/3  -  n]phltl 

1  '«(p2  2,  -  />,.,,  -  pi  ,,| 

Pl.iq  =  o  =  [/a  -  2(0  +  Re  V)\pi.iq 

3  in  (  Pl,  ,,  +  Pl.iq  -  Pl.iq  -  P,  2,  | 

+  2(0  +  Re  B)p,.,q  -  VqPl, 

Pl.iq  =  0  =  (/«  —  4/3 |p,. ,,  t-  i«(p,.2,,  -  p2.„| 

Pl.iq  =  0  =  [  - 1 ( 2 A  -  6)  -  2/31/5,,,, 

+  *«(Pt.I,  -  Pl.iq  I  -  <lqP  12 
Pl.iq  =  0  =  1(5/5.,  ,,,  +  in  |  Pl.iq  -  Pl.iq )  -  Vqpl.l 
+  2(0  -  Re  B)p,jq  +  2(0  +  Re  fl)/52  J, 
Pl.4,  =  0  =  [  -*(A  -  5)  -  (0  +  0) )  Pl.iq 

I  to  (  Pl.iq  —  Pl.iq  ~  Pl.iq]  ~  Vq  Pll 

Pl.iq  =  0  =  [/(A  +  6)  -  (0  +  D')\pi.iq 

+  '“Iam?  3-  A,.,,  -  Pl.iq]  -  VqPn  (3.13) 

where 

5  u iq  ro„,  A  =  u)  to,, 

and 

1,  =  *yVf’  (3.14) 

For  the  density  matrix  elements  p,„  3,  where  >>r  ^  3,  we 


gel  three  algebraic  conations: 

Pl.iq  =  0  =  f  -i(h  -  6)  -  3/3  +  B*]p,  3q 

~  'a Pl.iq  ~  rf?" 

Pl.iq  =  0  =  j/5  -  2(0  +  V.n  0)]p2.3,, 

-  in\Pu.q  I  Pl.iq  |  -  tqP'll 

Pl.iq  =  o  =  [»(A  +  6)/54.j,]  -  (0  -  B*)a,..v 

-  XqPi.x  -  in  Pl.iq  (315 

where 

iq  =  (3.  K 

For  the  density  matrix  elements  p  , ,,,  where  n  =  I,  •  •  • 
4,  we  get  four  algebraic  equations: 

Pi.iq  =  0  =  | '(A  I  fi)  -  3/3  t-  «]/>,.,,  +  ictpi.iq 

Pl.iq  ~  o  =  (<5  -  2(0  -  Ini  0)]/5,  j, 

3  to  [  p]  iq  4 

Pi.iq  ~  0  =  | iS  —  2(0  —  Re  0)]/5,  j,  =»  p,  ,,  =  0 

Pl.iq  =  0  =  (  “  f  ( A  -  5)  -  (0  -  0)]p,  4, 


3  inpi.lq  -  i'qPi.i 


(3.17) 


For  the  stationaiy  state,  we  get  />•,  3,  =  0  so  that  (3.15) 
arc  independent  of  (3. 17).  For  the  genera)  case,  we  have 
for  />„,,,,  15  algebraic  equations  which  can  be  solved  by 
simple  numerical  analysis  (excluding  p3  3,  =  0). 

The  use  of  the  density  matrix  elements  derived  in  the 
present  section  for  calculating  the  resonance  fluorescence 
spectra  is  described  in  the  next  section.  We  develop  an 
analytical  solution  for  the  resonance  fluorescence  spectra 
for  two,  two-level  atoms  which  are  at  resonance  (A  =  0) 
and  at  a  distance  smaller  than  a  wavelength. 

IV.  Resonance  Fluorescence  in  a  Strong  Field 

I  hc  resonance  lluorcsccncc  spectra  arc  given  as 
f  d 

**  ./  =  J'  I  Plq.lq  3"  Piq.iq  3  Plq.lq  ]  spoll,  (41) 

where  only  the  growth  due  to  the  spontaneous  emission 
(evaluated  by  Ihc  interaction  //’ )  can  be  seen  by  the  de¬ 
tector,  which  is  in  a  direction  different  from  that  of  the 
applied  field.  In  (4.1),  we  traced  over  the  atomic  states 
where  we  have  taken  into  account  that  there  is  no  spon¬ 
taneous  emission  to  the  upper  state  |</>,  >. 

Hy  using  the  equation 


=  -rlM'.P  1 

spool  ft 


we  get 


(/,  I  Plq.lq  !  sp»tni  =  ~2  Rc  IV>..2,I 


(4.3) 
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lit  l  I  sponl  ^  1  WqPlAil  I  ^  Re  l{,(P.l.  li)| 

(4.4) 

~h  [ ^v-^Kpom  —  — ^  IfijPi.vl'  ('*•-"’) 

Hie  four  terms  on  (lie  right  side  of  (4.2)-(4.4)  corre¬ 
spond  to  the  transitions  described  schematically  in  (lie 
diagram  of  Section  II.  In  order  to  obtain  the  resonance 
fluorescence  spectra,  we  need  explicit  expressions  for  (lie 
density  matrix  elements  p,„  ,,,,  which  can  he  calculated  by 
solving  independent  sets  of  algebraic  equations  given  in 
the  previous  section  as  (3. 10),  (3. 13),  (3. 15),  and  (3. 17). 
Since  the  treatment  of  the  general  cases  needs  the  use  of 
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simplifies  the  expressions.  The  condition  e  «  I  refers 
to  the  approximation  of  a  strong  field  where  the  Rabi  line 
width  0(7/2  is  large  relative  to  the  spontaneous  decay  line 
width  (i.  Tor  exact  resonance,  |5|  =  |  -  to  |  and  the 

approximation  |  5  | / /3  »  I  excludes  the  line  center  so 
that  our  approximations  are  valid  for  the  line  wings.  The 
solution  of  the  line  shape  under  these  approximations  by 
following  the  present  methods  is  interesting  as  it  can  be 
done  analytically  and  the  results  can  be  compared  to  those 
of  other  investigators  |7]-[13|. 

By  following  the  above  approximations,  we  have  ob¬ 
tained  solutions  to  p,„  ,,,  of  (3.13),  which  are  correct  to 
first  order  in  e.  By  substituting  these  solutions  in  (4.1), 
(4.3),  and  (4.4),  we  get  the  following  expression  for  the 
line  shape: 


(ST  -  2)  (6  I  4 SI ) c 


2S13(S23  -  8)  ( 3SI2  -t-  4)e 


[S2(SI2  -  2)1  +  c7( 8SI2  -  6)  1  SI(SI2  -  2 ) ( II2  -  8)]  +  e2(  I4SI4  -  82123  +  48) 


(4.9) 


very  lengthy  and  tedious  algcbia.  we  piclei  to  analyze  and 
discuss  these  cases  by  the  use  of  numerical  analysis  in  a 
separate  paper.  Here  we  limit  the  discussion  to  specific 
cases  for  which  the  theoretical  analysis  can  be  simplified 
and  for  which  the  physical  results  obtained  by  us  can  be 
compared  to  (hose  of  other  investigators. 

Wc  investigate  here  the  resonance  fluorescence  spectra 
for  a  very  stiong  external  field  at  exact  resonance  so  that 
instead  of  using  the  first  set  of  algebraic  equations  for 
p,„„  given  by  (3.10),  wc  use  the  approximate  solutions 
given  by  (3.11).  liven  under  (his  approximation,  the  an¬ 
alytical  solutions  of  p,„  „n  me  still  vciy  complicated  and 
tedious.  Therefore,  wc  limit  the  discussion  to  the  case  for 
which  kijj  «  I .  Tor  this  case,  H  =  (3,  e'(  '■  =  ca  ''  = 

I .  Assuming  that  the  two  dipoles  of  the  two  atoms  arc 
equal,  fj  =  0  and  (3.15)  and  (3.17)  arc  homogeneous 
linear  equations  whose  only  solutions  aic 

P  l.v»  ~  P2.  Irf  -  Pi.'r,  =  P. l.ly  ~  P  L’.|  “  P'.l'l  ~  1*  (d.6) 

1  hcicforc,  the  contributions  to  the  resonance  lltioieseencc 
spectra  from  the  density  matrix  elements  p ,  ,,  and  of 
(4.4)  and  (4.5)  vanish.  This  result  is  obvious  as  two  at¬ 
oms,  which  arc  within  a  distance  smaller  than  a  wave¬ 
length,  are  trapped  in  the  state  |  <A,  >  with  no  transitions 
to  other  states. 

Wc  introduce  the  definitions 

l>!  u  .  ‘"•I. 

or  or  n 


We  have  developed  an  analytical  expression  lor  the  line 
shape  for  the  special  case  whcic  the  following  approxi¬ 
mations  arc  valid: 


kr,i  «  I ;  A  =  0; 


*1 

ft 


»  I.  (4.8) 


Hie  assumption  of  exact  resonance  (A  =0)  and  the 
validity  of  the  relation  B  -  ft  for  kr{j  <■:.  I  very  much 


In  Tig.  I,  wc  rlcsciibe  the  line  shape  IT,  as  a  function  of 
the  normalized  frequency  52  =  |  oj,  —  |/or  for  various 
values  of  r  =  (1/ct.  Hie  line  shape  is  given  in  relative 
units  obtained  by  exact  calculations  of  the  expression  in 
the  curled  brackets.  The  maximum  oT  tbe  line  shape  at  52 
=  'll  corresponds  to  u,  —  co(l  =  a  -Jl  where  a  v/2  is  the 
Rabi  hequency.  This  maximum  becomes  narrower  for 
small  r  and  its  height  is  tending  to  °o  for  e  tending  to  zero. 
One  can  conclude  this  result  analytically  also  ns  the  sec¬ 
ond  term  of  (4.9)  becomes  singular  Tor  e  -*  0  at  12  =  'll. 
These  results  ate  not  surprising  as  they  are  in  agreement 
with  the  conclusions  made  by  other  investigators  and  the 
cflcct  is  similar  to  that  obtained  for  one  two-level  atom. 
On  the  other  hand,  our  conclusions  about  the  additional 
side  bands  are  different  from  those  of  other  authors.  Al¬ 
though  our  calculations  show  an  additional  maximum  near 
12  =  2  -Jl  (corresponding  to  twice  the  Rabi  frequency: 
2  -Jin  —  |  u,  -  co0  | ),  this  maximum  is  quite  broad  and  it 
docs  not  show  a  singular  behavior  for  e  tending  to  zero 
1 7 ) .  Taking  into  account  the  logarithmic  scale  used  in  Fig. 

I .  our  conclusion  is  that  the  additional  side  bands  are  quite 
weak  and  arc  not  in  agreement  with  the  relatively  strong 
resonances  predicted  by  other  authors  [7J-|  I0|,  1 13). 

Wc  can  explain  our  result  by  using  a  dressed  state  pic¬ 
ture.  The  additional  side  bands  result  from  transitions  in 
which  the  dipole  moment  is  doubled.  Hie  probability  for 
this  kind  of  transition  fora  stationary  two,  two-level  atomic 
system  is  quite  small,  as  the  relative  phases  of  different 
dipoles  lend  to  be  randomized.  In  a  recent  paper  1 12),  we 
have  suggested  observing  the  additional  side  bands  in  a 
time  dependent  spectra  where  a  strong  transient  coopera¬ 
tion  will  be  induced  by  initial  conditions. 

V.  Conclusions 

flic  theoretical  analysis  made  in  the  present  paper  for 
the  resonance  fluorescence  spectra  of  driven  two,  two-level 
atoms  can  be  used  for  evaluating  the  spectra  in  numerical 
calculations.  The  method  is  based  on  calculations  of  den- 
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Fig  I  The  line  shape  is  described  in  relative  units  accoitling  to  (4.9) 
as  a  function  of  the  normalized  frequency  II  ~  1 6 1 /or  =  J a*v  -  w, ,(/<». 
The  curses  denoted  by  (a),  (b),  and  (c)  correspond  to  normalized  spun 
taneous  broadening  *  =  0.1,  f  «  0.01,  and  t  -  0.001,  respectively. 
This  figure  represents  the  line  shape  under  the  approximations  given  by 
(4.8). 

sity  matrix  elements  p,„  „  and  p,„  m(  between  atomic  and 
atomic  plus  one  photon  stales.  We  have  presented  the 
equations  for  the  case  in  which  the  external  radiation  is 
propagating  perpendicular  to  the  interatomic  axis,  and 
have  introduced  the  explicit  spatial  dependence  of  the  co¬ 
herent  dipole-dipole  interaction  in  the  damping  constants 
for  the  density  matrix  elements.  The  equations  lor  p,„  „ 
and  p„,  nq  are  reduced  by  the  present  model  to  algebraic 
equations  given,  respectively,  by  (3.10)  for  p„,  „  and 
(3.13),  (3.15),  and  (3.17)  for  p,„  „9.  These  density  matrix 
elements  are  used  in  (4.1),  (4.3),  (4.4),  and  (4.5)  for 
evaluating  the  resonance  fluorescence  spectra.  For  the 
case  of  two,  two-level  atoms  at  resonance  interacting  with 
a  resonant  external  field,  which  are  within  a  volume  with 
a  dimension  smaller  than  a  wavelength,  we  have  derived 
an  analytical  expression  for  the  line  shape,  presented  in 
(4.9).  We  have  obtained  for  this  case  side  hands  shifted 
from  the  line  center  by  the  Rabi  frequency.  Additional 
side  bands  appear,  which  arc  shifted  ftom  the  line  center 
approximately  by  twice  the  Rabi  frequency.  The  addi¬ 
tional  side  bands  show  very  broad  and  weak  maxima,  and 
the  result  is  in  disagreement  with  the  singular  behavior 
predicted  in  previous  works  (7|-[l()|,  |I3|. 
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ABSTRACT 


The  quantum  fluctuations  in  intrinsic  bistability  of  a  two-level  system  are 
calculated  by  using  the  Einstein  relations.  All  the  second  moments  of  the 
dipole  operators  are  calculated  and  the  statistical  properties  of  the  reaction 
field  are  studied.  We  find  that  a  quantum  noise  which  is  above  the  standard 
quantum  limit  is  inherent  in  the  reaction  field. 
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1.  Introduction 

Intrinsic  optical  bistability  (IOB)  that  does  not  require  external  optical 
feedback  such  as  mirrors  has  been  described  in  various  theoretical  and  exper¬ 
imental  works.1-14  Most  of  the  observations  and  the  theoretical  discussions 
have  been  made  for  complicated  systems,  e.g.  in  semiconductors,  for  which 
a  full  quantitative  analysis  is  difficult.  A  relatively  simple  pure  quantum  me¬ 
chanical  system  which  produces  IOB  occurs  for  a  dense  collection  of  two-level 
atoms  interacting  via  the  electromagnetic  field  and  driven  by  an  external  ap¬ 
plied  coherent  field.  For  low  values  of  the  incident  field  the  strong  reaction 
field  reduces  the  net  field  into  a  low  transmitivity  branch.  As  the  incident 
field  is  increased  at  a  certain  critical  intensity  the  cooperative  interatomic 
interaction  is  broken  and  the  contribution  of  the  reaction  field  is  suddenly 
negated.  This  phenomenon  of  a  sudden  transition  from  the  low  transmitivity 
branch  into  the  high  transmitivity  branch  occurs  as  a  first  order  phase  tran¬ 
sition  far  from  thermodynamic  equilibrium.  Although  various  properties  of 
the  coexisting  steady  states  of  the  transmitted  radiation  have  been  described 
in  our  previous  works,  the  influence  of  this  phenomenon  on  photon  statistics 
has  not  yet  been  explored. 

The  purpose  of  the  present  work  is  to  study  the  quantum  fluctuations 
in  intrinsic  bistability  including  especially  the  statistical  properties  of  the 
reaction  field. 

In  a  recent  work15  we  have  treated  quantum  noise  fluctuations  in  IOB 
two-lvel  system.  However,  the  analysis  has  been  made  only  for  a  system 
with  a  rapid  dephasing  (1\  <<  T\)  and  only  the  affect  of  the  quantum  noise 
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on  the  bistable  states  has  been  discussed.  In  the  present  work  we  treat  a  pure 
quantum  mechanical  two-level  system  in  which  the  relaxation  time  2\  of  the 
inversion  of  population  is  the  spontaneous  decay  time  and  the  relaxation  time 
for  the  dipole  moment  is  given  by  T2  =  2Ti.  By  using  the  quantum  theory  of 
noise  developed  by  Lax16  we  calculate  all  the  second  moments  of  the  dipole 
operators  and  study  the  statistical  properties  of  the  quantum  noise  in  a  pure 
IOB  system. 

It  has  been  shown  by  other  authors  that  squeezing17-19  can  be  produced 
by  optical  bistability20  but  this  possibility  has  been  studied  only  under  the 
approximation  of  “a  good  cavity  limit”,  i.e.  for  bistability  in  high  Q  cavity21. 
For  such  cases  the  medium  relaxation  is  much  faster  than  that  of  the  field 
and  then  the  atomic  variabls  can  be  adiabatically  eliminated.  We  investigate 
here  the  problem  of  photon  statistics  for  the  extreme  case  of  a  “bad  cavity 
limit”21  where  we  preclude  any  reflecting  boundaries  and  where  the  relaxation 
of  the  field  in  the  medium  is  much  faster  than  the  relaxation  of  the  material. 
For  this  case  the  statistical  properties  of  the  reaction  field  are  related  to  the 
statistical  properties  of  the  dipole  operators. 

A  single  atom  driven  by  an  externally  applied  coherent  field  emits  sponta¬ 
neous  photons  that  are  “antibunched”  i.e.  the  intensity  correlation  function 
is  in  the  quantum  domain  </2l(0)  <  1  for  a  null  delay.22-23.  This  nonclassical 
property  of  the  light  observed  in  resonance  fluorescence  of  single  atoms  oc¬ 
curs  due  to  the  reduction  of  the  state  of  the  atom  into  its  ground  state  after 
the  emission  of  the  first  photon  so  that  it  cannot  immediately  emit  a  second 
photon.  The  antibunching  phenomena  vanish  for  the  radiation  emitted  in 
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resonance  fluorescence  of  many  two-level  atoms  due  to  destructive  interfer¬ 
ence  effects.  We  will  show  in  a  similar  way  that  the  quantum  noise  is  inherent 
in  the  IOB  system  and  that  squeezing17-10  cannot  be  obtained  in  this  system 
(the  bad  cavity  limit). 

In  a  recent  work  antibunching  behaviour  has  been  observed  in  the  reso¬ 
nance  fluorescence  light  emitted  by  a  many  atom  source.24  The  constructive 
phase  matching  conditions  in  that  work  were  obtained  by  using  a  four  wave 
mixing  scheme  in  which  the  intensity  correlation  function  was  measured  be¬ 
tween  two  detection  channels  of  spontaneous  radiation  emitted  in  opposite 
directions.  Although  this  system  is  basically  different  from  that  of  the  present 
work  the  methods  developed  here  can  be  used  also  for  treating  that  problem 
and  we  hope  to  do  it  somewhere  else. 

In  Sec.  2  we  describe  shortly  some  properties  of  Squeezed  States  (SS) 
and  explain  our  method  for  checking  the  existence  of  squeezing.  In  Sec. 
3  we  calculate  all  the  moments  of  the  dipole  operators  in  IOB  by  using 
Einstein  relations18,25.  The  statistical  properties  of  the  transmitted  radiation 
are  discussed  and  we  find  that  squeezing  is  not  obtained  in  the  present  system. 
In  Sec.  4  we  summarize  our  results  and  conclusions. 


2.  Squeezing  and  Anomalous  Correlation  Functions 


The  nonclassical  properties  of  quantum  light28  are  related  to  the  normal 
ordered  normalized  second  order  correlation  function  defined  for  stationary 


light  as 


9W(r) 


<  E(-\t)E(~\t  +  T)£?(+)(t  +  r)E^\t)  > 
<  E(~){t)E(+)(t)  >2 


(1) 


where  E^  and  E 1  1  are  the  positive-frequency  and  negative- frequency  parts 
of  the  elctromagnetic  field.  The  non-classical  properties  of  S.S.  are  usually 
obtained  by  mixing  the  L.O.  (coherent  state)  and  the  S.S.  in  interference 
experiments.27-28  It  is  generally  assumed  in  these  experiments  that  the  L.O. 
is  much  stronger  than  the  S.S.  so  that  the  denominator  of  Eq.  (1)  is  equal 
approximately  to  |a|4  where  |a|  is  the  absolute  value  of  the  coherent  state 
amplitude.  The  numerator  of  Eq.  (1)  can  be  developed  up  to  the  second 
order  in  the  S.S.  operators  and  the  result  is  given  by29-31 


<7(2)(r)  -  1  = 


<:  A[£?W(Q£?W(t  -f  r)j  :> 


(2) 


where  ::  denotes  normal  ordering  and  is  the  quadrature  amplitude  of 
the  S.S.:32 


E[e\t)  =  E\-  \t)ei9  +  E[+\t)e~i9  (3) 

6  is  the  phase  angle  of  the  coherent  state  which  is  changed  in  the  interference 
experiment  to  obtain  an  optimal  quantum  noise  reduction.  The  quadrature 
amplitudes  are  defined  in  the  interaction  picture  of  the  L.O.  so  that  the  rapid 
dependence  on  the  frequency  of  the  L.O.  is  eliminated.  The  function  g^2\r)  is 
dimensionless  and  since  the  normalization  is  arbitrary  it  is  chosen  so  that  |a|2 
is  the  flux  of  the  L.O.  photons  (sec-1)  and  consistently  <:  A[£l|8^(<)Ej^(<  + 
r)]  :>  is  given  also  in  terms  of  the  flux  of  S.S.  photons. 
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Non-classical  light  is  defined  by  the  antibunching  property  g^(r)  <  l.26 
Similar  expressions  to  Eq.  (2)  have  been  used  extensively  in  various  works 
for  analyzing  quantum  noise  reduction  under  the  Shot  Noise  Limit  (SNL) 
in  homodyne  and  heterodyne  experiments19,29-31  or  explaining  other  quan¬ 
tum  effects  in  correlation  measurements  between  two  photodetectors33.  In 
particular  the  Fourier  transform  of  Eq.  (2)  is  known  as  the  squeezing  spec¬ 
trum  and  various  theoretical  methods  have  been  developed  for  calculating 
this  spectrum.32,34-36 

The  anomalous  correlation  function  for  the  S.S.  is  defined  as35 

J++(r)  =<:  A +  r)]  :>=  |J++(r)|c^  (4) 

and  the  ordinary  correlation  function  as3B 

J_+(r)=<:A[4")(0^.+)(«  +  2')]->  (5) 

By  using  these  definitions  Eq.  (2)  can  be  expressed  as 

<7(2)(r)  -  1  =  cos2{6  +  */2)[2|J++(r)|  +  J.+{r)  +  J+.(r)]  (6) 

+sin2(0  +  <£/2)[-2]J++(r)|  +  J_+(r)  +  J+-(r)], 

where  [J+^.(r)]*  =  J _ (r).  It  will  be  shown  later  that  «/_+(r)  and  J+_(r)  are 

for  IOB  real  positive  functions  of  r.  Since  g^2\r)  —  1  expresses  the  amount 
of  quantum  noise  this  noise  is  minimum  for  sin2(0  -f-  <f>/2)  =  1  and  maximum 
for  cos2(6  +  <f>/2)  =  1  with  harmonic  dependence  on  the  L.O.  phase  8. 

In  the  present  work  we  are  interested  in  the  statistical  properties  of  the 
reaction  field  in  IOB.  It  is  sufficient  to  calculate  <7^(0)  in  order  to  see  if  anti¬ 
bunching  effect  (<7^(0)  <  1)  occurs  in  the  system  which  would,  in  principle, 
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lead  to  squeezing  spectra.  The  condition  p^(0)  <  1  gives  according  to  Eq. 

(6)  the  following  criterion  for  obtaining  squeezed  spectra 

2|J++(0)|>  J_+(0)  +  J+_(0)  (7) 

for  null  delay  times.  Squeezing  is  obtained  when  the  relation  (7)  is  fulfilled 
since  by  adjusting  the  phase  9  of  the  L.O.  we  get  the  non-classical  regime 
pl2)(0)  <  1  and  the  quantum  noise  is  reduced  below  the  SNL  in  the  homodyne 
or  heterodyne  or  other  correlation  experiments.  The  meaning  of  inequality 

(7)  is  that  for  obtaining  squeezing  the  anomalous  correlation  function  of 
creating  or  annihilating  two  photons  should  be  larger  than  that  of  the  usual 
correlation  function  of  creating  or  annihilating  one  photon. 

The  role  of  anomalous  correlation  functions  in  relation  to  squeezing  has 
been  treated  in  various  works.31,36,37  In  particular  the  constructive  anti¬ 
bunching  effects  described  in  Ref.  24  were  obtained  by  designing  the  ex¬ 
periments  so  that  constructive  interference  is  obtained  for  the  anomalous 
correlation  function  J++(r)  and  a  destructive  interference  is  produced  for 
the  ordinary  correlation  functions  J+_(t)  and  J_+(r). 

By  using  definitions  analogeous  to  Eqs.  (4-5)  and  the  inequality  (7)  the 
criterion  for  obtaining  squeezing  in  IOB  is 

2|  <  6<T+0£cr+0  >  |  >  <  6<T-06<r+o  >  +  <  6cr+o6<T-o  >  ■  (8) 

We  have  chosen  to  treat  in  the  present  work  a  pure  quantum  system  since 
as  is  well  known  the  insertion  of  additional  damping  mechanism  spoil  the 
effects  of  squeezing.  However,  we  will  show  that  even  in  the  pure  quantum 
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mechanical  system  a  quantum  noise  exists  which  is  above  the  Quantum  Noise 
Limit  (QNL).  We  calculate  in  the  second  section  all  the  second  moments  of 
the  dipole  operators  and  examine  the  statistical  properties  of  the  reaction 
field. 

3.  Statistical  Properties  of  the  Quantum  Noise  in  IOB 

We  consider  a  thin  sample  of  a  large  number  of  two-level  atoms  coupled 
to  each  other  only  via  the  electromagnetic  field  and  driven  externally  by  an 
applied  coherent  field.  In  our  previous  works  on  IOB  in  two-level  systems 
we  have  calculated  for  various  cases  the  expectation  values  per  unit  volume 
of  the  slowly  varying  dipole  moments  <  tr±0  >  and  for  the  inversion  of 
population  <  <*z  >  -12'14  In  the  calculations  for  the  expectation  values  the 
contribution  of  the  Langevin  noise  terms  vanish.  In  the  present  work  we 
assume  the  relations 


d'io  =  <  ^±0  >  -I  <hj±o;  &z  =<  <rz  >  -f  Scrz 


(9) 


and  linearize  the  equations  for  small  perturbations  of  the  system  from  the  sta¬ 
tionary  state.  By  following  the  derivations  given  in  our  previous  works12-14 


including  the  Langevin  noise  terms  we  get 


d(6a+ q) 
dr 


+ 


[-/?/ 2  4-  t(A  -  e  <  <rz  >j<$<7+o 

[wfl/2  -  ie  <  <r+o  >J££*  4-  /(-)  <  crz  >  \ 

~j~  =  -  unt><r+o  -  uJ*n6&- o  -  2/(+^cr+o)-  2 <  <r_0  > 

(10) 
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We  neglect  here  the  multiplications  of  with  8crz,8(r± 0-  The  present 

approximations  are  not  valid  when  the  quantum  fluctuations  are  large  or 
when  the  system  is  unstable. 

We  assume  a  fully  quantum  mechanical  model  where  (3  =  4jP|2fc3/3fe  is 
spontaneous  decay  time.  P  is  the  dipole  matrix  element  and  k  =  i u/c.  A  = 
u)  —  is  the  detuning,  where  u>  is  the  frequency  of  the  applied  field  E,  and 
u>o  is  the  atomic  resonance  frequency,  w/j  is  the  Rabi  rate  associated  with 
the  applied  coherent  field  and  e  is  the  frequency  renormalization  constant 
derived  for  a  dense  two-level  system  in  our  previous  work  as 

7tT  71,5c3 

e=~4^r 

The  frequency  renormalization  stems  from  coherent  dipole  dipole  interaction 
between  atoms  that  are  within  a  volume  of  a  cubic  wavelength.  The  free  field 
operators  of  Eqs.  (10)  obey  in  the  plane  wave  limit  the  equations: 

<  /'♦’W/'-’M  >=  W.t  - 

<  /MW/(tV)  >  =  <  /'-’M/'-V)  >  =  <  >=  0(L2) 

where  r  is  the  retarded  time.  These  operators  constitute  fluctuating  noise 
terms  which  stem  from  the  quantum  field  coupling  to  atoms. 

We  follow  in  our  theoretical  treatment  of  the  quantum  noise  in  IOB  the 
quantum  theory  of  noise  developed  by  Lax]6.  We  rewrite  our  fundamental 
Langevin  Eqs.  (10)  in  the  form: 

~  G^Sft)  +  t)  (13) 
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where  for  our  cases  a,G  and  F  are  three  dimensional  vectors  given  by: 


ai  —  8a+ 0;  a.2  —  ®3  —  8&z 


(14) 


Cl  =  [-/3/2  +  *(A  ~  c  <  crz  >)]<5<t+o  +  [wJj/2  -  ie  <  <r+0  >]<Sd-«, 

C2  =  Gf ;  G3  =  —/38crz  —  u)ji6ct+o  —  0  ,  (15) 

and 

Fi  =  f[-]  <  <rz  >;  /'2  -  /(f)  <  <rz  >;  F3  =  -2/0  <  <r+0  >  -2/0  <  ^_0  >  .(16) 
F^(t)  are  MarkofRan  operators  and  according  to  Eq.  (12)  obey  the  equations 
<  Fp{t)F„{u)  >=  2  <  Dul/(a,  t)  >  8(t  -  u )  (17) 

where 


<  On  >  — <  D 22  >  —  <  D 12  >  —  <  Hi3  >  — <  D32  >  —  0; 

2  <  023  >=  -2/1  <  <tz  ><  <7_0  >;  2  <  Z72i  >=  P  <crz  >2, 

2  <  ZI31  >=  -2/1  <  cr+0  ><  <7*  >;  2  <  f?33  >  =  4  <  cr^.o  X  <r_o  > 


One  should  notice  that  <  D ^  >  is  not  symmetric  as  FM  and  Fu  do  not 
commute. 

By  using  the  MarkofRan  property  of  the  operator  FM(t)  Lax  has  shown16 
that  tlie  equation  for  the  mean  motion  of  the  product  of  two  operators  a ^ 
and  a„  obeying  Eq.  (13)  is  given  by: 


dt 


2  <  >1  <  G^d,,  >  f  <  a^Gf,  >  . 


(19) 
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In  the  steady  state  when  the  second  order  operators  do  not  change  with 
time  we  obtain  the  standard  Einstein  relations 

—  2  <  D nv  >steady  state  =<-  Gnav  >»*  +  <  anGv  >,,  .  (20) 

The  order  in  these  equations  follows  from  the  quantum  nature  of  our 
system  and  in  this  respect  it  is  different  from  the  equations  for  the  moments 
in  a  Markoffian  system  which  can  be  derived  classically. 

We  introduce  the  following  definitions: 

p  =  -0/2  +  i[A  -  e  <  at  >];  q  =  u)*n/2  -  ie  <  <r+0  >  .  (21) 

Einstein  equations  for  our  system  are  obtained  from  Eqs.  (20)  by  substi¬ 
tuting  the  values  of  a ^  Gv  and  DMU  according  to  Eqs.  (14),  (15)  and  (18), 
respectively,  and  by  using  the  short  notation  of  Eq.  (21). 

We  get  the  following  9  Einstein  equations  for  our  system: 

<  6az6a+0  >  +  <  Sv+oSa,  >=  -(2 p/q)  <  Sa+0Sa+0  >; 


<  Saz6a_ o  >  -I-  <  6cr_0Scrz  >-  -(2 p* jq')  <  £<r_0£cr_o  >; 


u>rp 


<  Sa\  >=  2  <  a+ o  ><  <r_0  >  <  6a+0Sa+0  > 
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+  <  6a-06a-o  >  ; 

0q 


q  <  6az6a_o  >  +9*  <  Sa+oSaz  >=  0  <  Sa+o^-o  >1 
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(p-/3)  <  8ct+o8<tz  >  +<j  <  8crz8oz  >=  u>n  <  6<t+o5(j+o  > 

<  <J«T+0<5cr_o  >  ; 

(p*  -  (3)  <  8c t-0S<tz  >  4 -q*  <  8<tz8ctz  >=  u i'R  <  «?<T_o^cr_o  > 

-\-u>r  <  >  +2/3  <  <rz  ><  <j_ o  >  ; 

q  <  8cr_08<Tz  >  +9*  <  8crz8o+0  >=  — /3  <  erf  >  +/3  <  <S<r_06er+0  >; 

(p  -  / 3)  <  <5<t*<5<7+0  >  +g  <  <5 erf  >  -wfl  <  6<r+06<r+0  > 

-wjj  <  <5<7_06<r+0  >=  2/3  <  tr+0  ><  <rz  >  ; 

(P*  -  0)  <  8crz8cr_o  >  +g*  <  6trf  >  -u>n  <  5<r+05er_0  > 

— <  8 <T _q8 o'  _ q  >=  0  (22) 

Eqs.  (22)  consist  of  9  algebraic  equations  from  which  one  can  calculate 
numerically  all  the  moments  for  our  system.  A  careful  analysis  of  these 
equations  shows  that  they  fulfill  the  following  relations: 

<  8oz8<j 4.0  >  =<  6<r_o<5cr,  >;  <  8ctz6(T- o  >*  =  <  ^er+0^erz  > 

<  6o+08o+0  >  =<  (5<r_o5cr_o  >  (23) 

The  functions  <  8o\  >,  <  <5cr+o6cr_0  >  and  <  6cr_05(r+ o  >  are  real  and 
positive  but  <  6<r+06er_0  >7^<  5er_06(r+o  >  due  to  the  unsymmetric  property 
of  D 
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By  using  the  above  relations,  Eqs.  (22)  can  be  reduced  to  the  following 
form 


2  p 

<  {5crz,5cr+o}  >= - <  6(r+08(r+o  > 

9 


(24~<t) 


<  8a]  >=2  <a+0><  a_0  >  <  8a+05a+0  >  +c.c.);  [21 -b) 


(3q 


( q *  <  {8az,8a+ 0}  >  +c.c.)  =  -<az>7/3  +  (3<  {8a+0,  8a„o}>;  (2.4  -C) 


(q,{8<rzy8a+ 0]  +  c.c.)  =  -0  <  crz  >2  +/?  <  [50-_o,<5<r+o]>;  (2-4'  d) 


(p  -  (3)  <  {6az,8a+0}  >  +2 q  <  8a]  >  = 

2u>n  <  6c+o8a+o  >  +wr{^<7’+0) ^^-o}  +  2/?  <  <x+0  ><  <7Z  >5 

(p  -  (3)  <  {8<jz,8a+ 0]  >=  u^[£<r_o,  ^+o]  +  2/0  <  <r+0  ><  cz  >; 

where  in  these  equations  {  }  and  j  ]  denote  anticommutation  and  commu¬ 
tation  relations  respectively. 

In  the  present  work  we  are  treating  the  fluctuations  in  the  reaction  field 
and  check  the  possibility  for  getting  squeezed  light  according  to  Eq.  (8). 
For  this  purpose  it  is  enough  to  calculate  the  moments  <  6cr+06(r+o  >=< 
8a-o8a_0  >*  and  <  6er+05cr_o  >  +  <  <5<r_06cr+0  >=  {6(r+0, 6cr_0}.  By  substi¬ 
tuting  Eq.  (24a)  into  Eq.  (24c)  we  get: 

<  {8a+,8a_0}  >  +~~  <  8aJ,08ay0  >  +2~  <  8a_06a_0  >  =  <  az  >2  (25) 

(in  /irt* 


(ll-e) 

(21-f) 
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By  substituting  Eqs.  (24-a)  and  (24-b)  into  Eq.  (24-e)  we  get 


<  {6(t+o,8o’_q 


>  +  < 


r--n 

-  <  ><:  >  -  2/3  <°>°"<  it  ft6) 

q*(3  4  4 


We  find  that  the  moments  needed  for  calculating  the  properties  of  the 
reaction  field  can  be  calculated  from  Eqs.  (25)  and  (26)  without  the  need  to 
evaluate  all  other  moments.  Let  us  introduce  the  following  definitions 


<  8<t+06<x+0  >  —  =  A\  <  {£<r+0,6<7--o}  >=  2 r; 
ur 


(27) 


where  r  is  real  and  positive  and  the  condition  for  squeezing  becomes  r  <  |j4|. 
For  the  steady  state  solutions  of  IOB  the  following  relations  are  fulfilled 

~(wr/2)  <(TZ> 


<  <r+0  >=  - 


<  <r+ o  ><  <r_ o  >  = 


*( A  -  e  <  <rr  >)  -  (/3/2)  ’ 
-  c  <rt  >  [<  az  >  +lj 


(28) 


Using  these  relations  the  real  and  imaginary  parts  of  Eq.  (26),  and  Eq.  (25) 
consist  of  three  algebraic  equations  from  which  the  parameters  r  and  A  can 
be  calculated  and  the  condition  r  <  [A|  for  squeezing  can  be  examined  by  the 
numerical  calculations.  The  phase  of  wR/uR,  fixed  by  an  arbitrary  phase  in 
the  external  field,  will  afTect  the  phase  of  <  6cr+06(r+0  >  but  will  not  change 
the  results  for  the  values  of  r  and  |v4|  =  |  <  6<r+06cr+0  >  |. 

By  substituting  the  value  of  <  cr+0  >  from  Eqs.  (28)  into  Eq.  (21)  we 

get: 


q  =  (<4/2)1 


(jA  -  (3/2) 

t(A  -  e  <  <rz  >)  -  /3/2 


(29) 
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Since  general  analytical  solutions  of  Eqs.  (25-26)  are  still  quite  compli¬ 
cated  we  pursue  further  our  quantitative  analysis  for  three  special  cases  and 
obtain  qualitatively  general  conclusions. 

Let  us  develop  the  analytical  solutions  under  the  resonant  condition  A  = 
0  where  we  get: 


q  =  (u>n/ 2)---  v  =  -0/2  -  te  <  at  >  (30) 

p  +  zie  <  a z  > 

Substituting  Eqs.  (27),  (28)  and  (30)  into  Eqs.  (25-26)  we  get 

2r  _  _  A-W  -  2“  <  >)*  =< >»  (3i) 

0(0  -  2 te  <  crz  >)  0  +  2ie  <  az  > 


„  ,  ,  0  +  2te  <  <rz  >  (0  +  2te  <  >)2(3/3  +  2ie  <  crt  >) 

2r  +  .4(2+ - - - + - ’’ 

+A'mT  2V:  <  ~  =  a — oT  “ - 13  <  v.  >!  +  <  v,  >|.  (32) 

0(0  -f-  2ie  <C  crz  >)  /9  4-  2tc  <  crz  > 

We  study  the  quantum  fluctuations  in  the  reaction  field  both  for  the  ordinary 
Bloch  equations  (e  =  0)  and  for  the  IOB  case  (e  >>  0). 

For  the  ordinary  Bloch  equations  on  resonance  we  have 

e  =  0;  02 /\u>n\2  =  -2  <  <rz  >  /[<  <rz  >  +1]  (33) 

and  Eqs.  (31-32)  can  be  written  as 

2r  —  (.4  +  /T)  =<  <rt  >2,  (34) 


2 r  +  3,4(1  +  p— -)  +  A*  =  -3  <  <rt  >2  -  <  <rt  > 

I 
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After  a  straightforward  algebra  we  get 


Aj  =  ImA  =  0;  Ar  =  Re  A  =  <  <rz  >2  [<  <rz  >  +1]  — 

3  6 

r  =  ReA+  <  az  > 2  /2.  (36) 

In  Fig.  1  we  describe  the  values  of  r  and  |A|  as  a  function  of  <  crt  > 

according  to  Eqs.  (36)  for  Bloch  equations  on  resonance.  We  find  that  r  >  |A| 

so  that  the  reaction  field  is  bunched  and  squeezing  cannot  be  obtained  for  the 

transmitted  radiation.  We  find  also  the  conclusion  that  the  fluctuations  in 

the  ordinary  correlation  function  are  decreasing  as  a  function  of  the  intensity 

of  the  external  field  so  that  they  are  negligible  in  the  region  of  <  <rz  >~  0 

and  maximal  in  the  region  <  <xz  >~  —1. 

Let  us  now  discuss  the  solutions  of  Eqs.  (32-33)  for  resonance  (A  =  0)  in 

intrinsic  bistability  (e  >>  (3).  For  this  case  the  relation 

|  it  _  [<  <rz  >  +l)[/?2  +  4e2  <  <rz  >2} 

|W/l1  2  <  <rz  > 

is  fulfilled.  We  simplify  further  Eqs.  (32-33)  under  the  approximation 

2e|  <  <rz  >  |  >>  (3  .  (38) 

This  approximation  is  very  good  for  the  lower  branch  of  intrinsic  bistability 
where  we  have  relatively  stronger  quantum  fluctuations.  We  get  under  this 
approximation: 

2r  +  6AR  -  Ae<°‘>Al  =<  <r,  >2  (39) 


2r  T  Ar 


10  <  az  > 

<  crt  >  +1 


a  4e  5  >2 

1  (3\<  <TZ  >  -f  1| 


4e2  <  <tz  > 2 


[3  <  <rz  > 2  +crt  >j(40) 
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(41) 


.  .  .  10  <  <rz  >  4e  <  <rz  >  4e  <  <rx  >2 

+  +  /,»l— +  j3[<  <T.  >  +11 


2e  <  crz  > 


(3  <  trx  >2  +  <  crt  >] 


In  Fig.  2  the  values  of  r  and  j.4)  are  calculated,  as  a  function  of  |  < 
crz  >  |  >  0.5  according  to  Eqs.  (39-41)  for  an  IOB  example  on  resonance 
(A  =  0)  for  which  e  =  10 0  and  the  condition  (38)  is  valid.  We  find  that  the 
quantum  fluctuations  for  an  IOB  system  are  similar  to  those  of  the  ordinary 
Bloch  equations.  Similar  to  the  ordinary  Bloch  equations  r  >  j>4 j  so  that 
squeezing  is  not  obtained.  For  the  upper  branch  in  IOB  the  terms  which 
are  proportional  to  e  <  <rz  >  are  quite  small  so  that  the  equations  become 
similar  to  the  ordinary  Bloch  equations  where  squeezing  is  not  obtained. 

Let  us  now  discuss  as  a  third  case  the  solutions  of  Eqs.  (25-26)  for  IOB 
under  the  resonant  condition  A  —  e<<rt>=0  where  the  relation 

02/Wn\2  —  _2  <  <rz  >  /[<  <rz  >  +1)  (42) 

is  valid  (as  in  Eq.  (33)).  For  this  case  the  relations 

p  =  -0/2]  <  <r+0  >= - - - ;  q  =  (wr/2)[- - - - )  (43) 

are  fulfilled.  By  using  Eqs.  (43)  and  Eqs.  (27)  we  can  write  Eqs.  (25-26)  as 

2r - 1  V'2.< < »;> + c-c-i”<<r->  -  <44> 


2r  +  M|3+;ML-_J - ]  +  A-0 

\u>n\2  0  -2ie  <  ( rt  >  0  -V  2te  <  <rz  > 
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Using  the  relation  (42)  and  the  approximation  (38)  we  reduce  Eqs.  (44-45) 
into  the  form: 


2  r  +  2Ar  4-  Aj- 


0 


e  <  or,  > 


=  <  O’,  >“ 


(46) 


2 r  +  Ar[ 2  - 

=  -2  <  crz  > 


[<  <rz  >  +  l]2e2  <  <rz  +  ^ 

2  -  <  <TZ  >  [<  <TZ  >  + 1]  , 


3/3 

e[<  az  >  +1) 


0 


e  <  crz  > 


(47) 


AM 


3  (32 


[<  az  >  +l]2e2  <  <tz> 


+ 


3/3 


e  <  az  >  [<  crz  >  +l]e 


2  <  <rz  >2  e[<  <rz  >  +1] 

0 


(48) 


In  Fig.  3  the  values  of  r  and  \A\  are  calculated  as  a  function  of  |  < 
(tz  >  |  >  0.5  according  to  Eqs.  (46-48)  for  an  IOB  example  on  the  resonance 
A  —  e  <  az  >  =  0  for  which  e  =  10/3  and  the  condition  (38)  is  valid.  We  find 
that  the  quantum  fluctuations  in  this  case  are  very  large  relative  to  those  in 
Fig.  2.  Squeezing  is  not  obtained  in  this  example  and  the  reaction  field  is 
again  bunched.  Although  we  evaluated  the  parameters  A  and  r  according 
to  Eqs.  (25-26)  under  special  conditions  we  can  obtain  general  conclusions. 
Under  the  resonance  condition  A  —  e<<7z>=0  the  absolute  value  of  the 
parameter  in  Eq.  (25)  is  minimal.  On  deviation  from  this  resonance 
we  would  expect  therefore  a  decrease  in  the  parameter  |<  5<r+0£<r+0  >|,  i.e. 
worse  conditions  for  squeezing,  We  summarize  our  conclusions  in  the  next 
section. 


4.  Results  and  Conclusions 


The  statistical  properties  of  the  reaction  field  in  IOB  has  been  studied. 
The  quantum  properties  of  the  reaction  field  are  related  to  the  anomalous 
correlation  function  and  a  general  criterion  for  obtaining  squeezing  is  given 
(Eq.  (8)]. 

General  equations  for  calculating  all  the  second  moments  of  the  dipole 
operators  have  been  developed  by  linearizing  the  equations  for  small  per¬ 
turbations  of  the  system  from  the  stationary  state  and  by  using  the  corre¬ 
sponding  Einstein  relations  [Eqs.  (22-24)].  The  fluctuations  in  the  reaction 
field  can  be  calculated  numerically  from  Eqs.  (25-27)  where  the  condition 
2  |<  £<r+0(5cr+0  >|>  <  8a+08a_0  >  +  <  8<t_08(t+0  >  represent  the  antibunch¬ 
ing  property  (and  vice  versa). 

By  solving  our  equations  for  special  examples  and  examining  the  general 
behaviour  of  the  equations  we  find  that  the  reaction  field  is  bunched  so  that 
a  quantum  noise  which  is  above  the  quantum  noise  limit  is  inherent  in  the 
IOB  system.  We  may  conclude  therefore  that  the  photon  statistics  in  the 
bad  cavity  limit  is  more  chaotic  than  that  of  the  good  cavity  limit  in  which 
squeezing  was  predicted.  We  find  however  the  conclusion  that  the  quantum 
fluctuations  in  the  steady  states  of  IOB  on  resonance  (A  =  0)  are  similar  to 
those  predicted  for  the  ordinary  bloch  equations  (e  =  0)  on  resonance  for  the 
same  state  of  inversion  of  population. 

The  calculations  made  in  this  work  show  that  the  quantum  noise  in  IOB 
of  a  pure  quantum  mechanical  system  is  quite  large.  By  inserting  large 
dephasing  processes16  we  can  improve  the  conditions  for  IOB  but  not  for 
obtaining  squeezing. 
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FIGURE  CAPTIONS 


Fig.  1.  The  parameters  r  and  |A|  are  described  as  a  function  of 

<  crz  >  for  the  ordinary  Bloch  equations  on  resonance  (A  =  0).  The 
calcualtions  have  been  made  according  to  Eqs.  (36). 

Fig.  2.  The  parameters  r  and  \A\  are  described  as  a  function  of  az 
(for  | az\  >  0.5)  for  IOB  on  resonance  (A  =  0)  where  e  =  10/3. 

Fig.  3.  The  parameters  r  and  |A|  are  described  as  a  function  of 

<  az  >  (for  |<r2|  >  0.5)  for  IOB  in  which  A  —  e  <  crz  >—  0  and  e  =  10/3. 
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